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ABSTRACT 

'  An  efficient  method  is  presented  for  calculating  equilibrium  chemical  composition 
at-  a  given  pressure  and  stoichiometry  with  temperature,  enthalpy,  or  entropy  specified 
as  the  additional  condition.  A  substantia]  reduction  in  computational  time  was  achieved 
over  conventional  methods  by  employing  a  linearization  technique  and  reducing  the  number 
of  equations  to  be  solved  to  a  small  number  of  linear  algebraic  equations.  Examples  are 
given  for  the  H2/O2  system  (six  equations  reduced  to  two)  and  C/H2/02/N2/Ar  system 
(eleven  equations  reduced  to  three).  A  computer  program  is  presented  for  solving 
problems  with  fuels  which  involve  C  and  H  (hydrocarbon  or  hydrogen)  or  N  and  H 
(such  as  hydrazine)  and  with  oxidants  such  as  O2,  H2O2,  and  HNO3.  Extension  of  the 
method  to  include  other  systems  is  explained  in  sufficient  detail  to  eliminate  the 
requirement  for  an  extensive  knowledge  of  chemistry  by  potential  users. 
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NOMENCLATURE 

aji  Matrix  coefficient 

a,b,c,x,y,z  Relative  atom  moles  of  nitrogen  (N),  argon  (Ar),  carbon  dioxide  (CO2), 
carbon  (C),  hydrogen  (H),  and  oxygen  (O) 

bik  kth  coefficient  in  polynomial  thermodynamic  property  of  species  i 

Cp  Specific  heat  of  mixture,  cal/°K 

Cp  Mean  specific  heat,  Btu/lb°R 

CPi  Specific  heat  of  species  i,  cal/mole°K 

dji  Matrix  coefficient 

F  Gibbs  free  energy  of  the  mixture,  cal/mole°K 

Fx.i  Standard  Gibbs  free  energy  of  species  i  at  1  atm  and  temperature  T, 

cal/mole°K 

H  Enthalpy  of  the  mixture,  cal/gm 

H  Mean  enthalpy,  Btu/lb 
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hj  Enthalpy  of  species  i,  cal/mole°K 

J(ajj)  Determinant  of  matrix  of  coefficients 

K  Equilibrium  constant 

m  Number  of  independent  species  which  determines  order  of  matrix 

n  Total  number  of  species 

P  Pressure,  atm 

Pi  Partial  pressure,  atm 

Pv  Vapor  pressure,  atm 

Qcomb  Heat  of  combustion  of  a  fuel,  Btu/lb 

q  Number  of  elements 

R  Gas  constant,  cal/mole°K  or  cm3  atm/mole°K 

r  Ratio  of  number  of  atom  moles  of  a  species  to  that  of  another 

S  Entropy  of  mixture,  cal/gm°K 

Sf.i  Standard  entropy  of  species  i  at  temperature  T  and  1  atm,  cal/mole°K 

T  Temperature,  °K 

XX  Dummy  variables 

5  Decrement  of  specie  i 

Stoichiometric  coefficients  of  species  i  in  reaction  j 
<j>  Equivalence  ratio 

SUPERSCRIPTS 
o  Standard  state 

*  Specified  input  value 

'  Current  estimate  of  variable 
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SUBSCRIPTS 

ij,k,/3  Dummy  indices 

i  Property  of  ith  species 

j  Property  of  jth  reaction 

ref  Property  evaluated  at  reference  temperature 

T  Temperature 
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SECTION  I 
INTRODUCTION 

There  are  many  problems  involving  the  performance  of  propulsion  systems  in  which 
it  is  necessary  to  know  the  equilibrium  properties  of  combustion  gases.  As  long  as 
combustion  temperatures  are  relatively  low  (T  <  3000  to  3500°R)  good  results  can  be 
obtained  by  assuming  no  dissociation  of  the  combustion  products  and  including  an 
empirical  correction  for  dissociation  as  a  function  of  temperature  only.  An  effective 
isentropic  exponent  can  be  determined  and  fluid  dynamic  properties  may  be  calculated 
using  conventional  aerodynamic  functions. 

When  it  is  necessary  to  evaluate  systems  where  higher  temperatures  occur,  such  as 
turbojets  with  afterburners,  ram  rockets,  or  ramjets,  low  temperature  empirical  corrections 
rapidly  become  very  poor.  A  much  more  difficult  equilibrium  chemistry  calculation  must 
then  be  carried  out  since  the  complete  set  of  nonlinear  equations  which  describe  chemical 
equilibrium  must  be  solved. 

Solution  of  the  complete  set  of  nonlinear  equations  is  very  laborious  by  conventional 
methods,  even  for  a  relatively  simple  system  such  as  a  laboratory  flame.  Computers 
alleviated  this  problem  (Refs.  1  through  3).  Increasing  system  complexity  and  requirements 
for  thousands  of  data  point  calculations  brought  new  problems  directly  associated  with 
the  old,  time-consuming  methods  of  calculation. 

Numerical  techniques  for  solving  systems  of  nonlinear  equations  have  been  developed 
to  a  high  degree;  therefore,  it  is  necessary  to  devise  operations  in  the  system  of  equations 
to  effect  significant  reductions  in  the  overall  calculation  time.  The  most  successful  operation 
which  has  been  applied  is  that  of  linearization,  using  a  truncated  Taylor  series  expansion 
about  the  solution  (Refs.  4  and  5).  The  method  of  Ref.  4  employs  an  expansion  of  the 
free  energy  of  the  mixture,  using  the  method  of  steepest  descent  to  minimize  this  quantity 
corresponding  to  chemical  equilibrium.  The  method  of  Ref.  5  employs  an  expansion  of 
the  conservation  equations  and  the  Newton-Raphson  method  to  solve  for  the  equilibrium 
composition. 

The  technique  described  in  this  report  is  essentially  that  of  Weinberg  (Ref.  5)  extended 
to  allow  temperature  to  be  calculated  simultaneously  with  the  chemical  composition  and, 
in  addition,  to  include  a  numerically  stable  computational  technique.  When  temperature 
is  to  be  determined,  either  enthalpy  or  entropy  may  be  specified  at  a  given  pressure, 
thus  allowing  the  temperature  to  be  obtained  from  an  energy  balance  on  an  engine  or 
from  the  known  entropy  at  some  point  along  an  isentropic  expansion  or  compression. 

SECTION  II 

CALCULATION  OF  EQUILIBRIUM  COMPOSITION 

In  order  to  calculate  the  equilibrium  composition  of  a  mixture  of  n  components 
(species)  at  a  given  temperature,  pressure,  and  equivalence  ratio,  one  needs  n  independent 
equations.  For  a  system  containing  q  elements  from  which  n  species  are  constituted,  it 
is  possible  to  write  q  conservation  equations  which  are  linear  combinations  of  the 
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concentrations  of  species,  and  n-q  chemical  equilibrium  equations  which  are,  in  general, 
nonlinear  combinations  of  the  species  concentrations.  These  chemical  equilibrium  equations 
define  the  relative  amounts  of  each  species  at  a  particular  pressure  and  temperature. 

In  general,  it  is  possible  to  choose  any  q  species  which,  as  a  group,  contain  all  the 
elements  and  write  equilibrium  equations  for  the  remaining  species.  Substituting  these 
equations  into  the  conservation  equations  results  in  a  set  of  q  nonlinear  equations  in  q 
unknowns. 

Several  methods  have  been  developed  for  solving  nonlinear  equations  with  the 
Newton-Raphson  method— the  most  frequently  used.  There  are  limitations  to  the  use 
of  this  method,  but  it  is  a  powerful  tool,  and  convergence  can  usually  be  guaranteed. 
The  nonlinear  equations  are  linearized  by  expanding  the  equations  in  a  first-order  Taylor 
series.  In  order  to  solve  the  linearized  set  it  is  necessary  to  have  approximate  values  of 
concentration  for  q  of  the  species.  The  set  of  q  linearized  equations  is  then  solved  for 
corrections  to  the  approximate  values.  The  process  is  repeated  to  find  new  corrections 
to  the  previous  values.  Although  the  Newton-Raphson  technique  always  converges  to  a 
solution,  some  experience  is  necessary  to  ensure  that  the  physically  realistic  solution  is 
obtained.  In  addition,  the  rate  of  convergence  is  greatly  improved  if  experience  is  used 
to  specify  which  of  the  q  species  should  be  chosen  and  good  approximations  are  made 
for  the  initial  values.  Finally,  if  the  number  of  equations  can  be  reduced,  considerable 
savings  in  computational  time  are  obtained  since,  if  there  are  q  equations  for  q  unknowns, 
q2  operations  are  necessary  to  solve  for  the  q  unknowns. 

The  method  developed  in  this  report  is  essentially  the  application  of  experience  to 
improve  established  methods  for  equilibrium  calculations.  The  method  is  extremely  efficient 
in  obtaining  solutions  for  a  variety  of  complex  systems  including  the  simultaneous  solution 
of  the  energy  or  entropy  equations  to  determine  temperature  if  it  is  an  additional  unknown. 

The  species  are  assigned  partial  pressures  P;  and  the  necessary  equations  are 
(temperature  is  not  a  variable  in  the  following  derivation) 

Pi  =  g(P/g)  *  =  l,n-m;/3  =  1,  m  (1) 

where  the  functions  g(P)  are  the  equilibrium  equations  and  may  include  one  or  more 
of  the  conservation  equations.  The  conservation  equations  are  expressed  in  the  form 

ap  =  s(pK^  K  =  1»n5  Is  =  1> m  (2) 

where  a.p  is  defined  in  such  a  way  that  it  is  identically  equal  to  zero  when  the  solution 
has  been  obtained.  Expanding  Eq.  (2)  in  a  Taylor  series  and  neglecting  orders  higher  than 
the  first 


a/3  =  a/S 


n 

+  L 

K=l 


da 


/S 


dP, 


/8  =  1,  m 


(3) 
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Differentiating  Eq.  (1)  one  obtains: 


5.  mf  ^ 


j3=i 


«?p 


/s 


(4) 


Note  that  since  the  equations  are  to  be  solved  in  finite  difference  form,  finite  difference 
notation  is  used  throughout. 

Substituting  Eqs.  (4)  in  Eqs.  (3),  evaluating  the  partial  derivatives  using  Eqs.  (1), 
and  letting  a p  approach  0,  one  obtains  m  equations  for  the  m  unknown  5P.  The  Sp's 
are  the  desired  corrections  to  the  initial  values  of  Pp  which  are  denoted  by  Pp.  Thus 
the  corrected  value  is  given  by 

P/3  =  p/3  +  5/S  (5) 

and  the  remaining  (n  -  m)Pj’s  are  obtained  by  substituting  the  Pp's  in  Eqs.  (1).  The  process  is 
repeated  until  the  6p's  are  within  acceptable  limits. 

Some  examples  of  increasing  complexity  are  given  which  illustrate  this  method 
of  calculation  and,  in  addition,  a  method  is  included  which  simultaneously  solves  for  the 
temperature  when  enthalpy  or  entropy  is  known. 

2.1  APPLICATIONS  OF  THE  METHOD 


2.1.1  Hydrogen/Oxygen  (m  =  2) 

The  following  n  equations  completely  describe  the  relationship  between  the  n 
unknown  partial  pressures. 


Dissociation  Equilibria 


Species 


j  =  1 

0.5  H2  *7  H 

P,  -  Kj-v/Pj 

j  =  2 

0.5  02  -  0 

^2  =  K2\/P6 

j  =  3 

1120  *7  0.5  1I2  +  OH 

P3  =  K3P4A/P7 

j  =  4 

H20  H2  +  0.5  02  ■ 

p4  -  k4p5>/p7 

i  =  1  2  3  4 

5  6 

H  O  Oil  H20 

h2  o2 

Element  Conservation  Equations 


H  atoms: 

2PS  + 

P1  + 

P3  + 

2P4  = 

PH 

O  atoms: 

2Pg  + 

P2  + 

P3  + 

P4  " 

po 

(6) 


(7) 
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In  addition,  there  is  the  following  conservation  equation  for  species 

n 

I  pi  -  p  (8) 

i=l 


Let  the  ratio  of  H  atoms  to  O  atoms  in  the  unbumed  mixture  be  denoted  rn  o  =  pH  / po  ■ 
The  O  atom  conservation  equation  can  then  be  substituted  into  the  H  atom  conservation 
equation  as  follows: 

2P5  •*-  Pj  +  P3  +  2P4  -  rHO(2P6  +  P2  +  P3  +  P4)  =  0  (9) 

One  thus  has  six  equations  for  the  six  unknown  partial  pressures.  The  equilibrium  constants 
Kj  are  known  functions  of  temperature  (see  Appendix  I).  It  remains  to  choose  the  two 
species  to  which  the  remaining  four  species  are  to  be  related.  The  two  appearing  in  the 
largest  concentration  should  be  selected  in  order  to  minimize  numerical  problems.  At  low 
temperatures,  dissociation  will  be  almost  zero.  The  partial  pressure  of  atomic  hydrogen 
and  oxygen  and  hydroxyl  radicals  will  be  negligible.  For  fuel  lean  calculations  the  partial 
pressure  of  hydrogen  molecules  will  be  negligible  as  will  be  the  partial  pressure  of  oxygen 
molecules  for  fuel  rich  calculations.  For  stoichiometric  calculations,  all  the  partial  pressures 
are  small  except  that  of  water  vapor.  Thus,  water  vapor  is  a  natural  choice  and  hydrogen 
and  oxygen  are  used  for  rich  and  lean  calculations,  respectively.  Stoichiometric  calculations 
are  made  with  hydrogen  as  the  second  variable  although  oxygen  could  be  used  instead. 


FUEL  LEAN  CASE 

-V  V  .»/ 

P,  = 

P2  =  K2P62  -. 

These  equations  represent  those  expressed  as 


p3  .  K3kHP6 

P5  ‘ 

Eqs.  (1).  Differentiating  these 


(10) 

equations 


Sj  =  0.5  ^-54  -0.25-^56 

4  6 

*2  =  °-5~TS6 
6 

§3  =  0.5  a4  +  0.25  ^  s6 

4  6 


s 


5  = 


(ID 


These  equations  represent  those  expressed  as  Eqs.  (4). 
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Let 


and 


a\  -  2P5  +  Pj  +  P3  +  2P4  +  rHO(2P6  +  P2  +  P3  +  P4) 


«» -  E  Pi  -  ? 

j=l 


(12) 


These  equations  represent  those  expressed  as  Eqs.  (2).  Differentiating  these  equations 


da  1 
*  1 

da  j 

dal 

5F,*1 

5p7  =  ~  rH0 

1 

J*1 

II 

da\ 

dal 

da  1 

rHO  d?l=  2 

(9a  9 

i  =  1,6 

rHO 


rHO 


(13) 


These  expressions  together  with  Eqs.  (11)  above  represent  those  expressed  as  Eqs.  (3). 

METHOD  OF  SOLUTION 

r 

1.  Assume  values  for  P4'  and  Pg'  using  the  stoichiometric  equation 

P4  -  20P/(1  +  0)  Pfi  =  P  -  P4  rH0  =  2<f> 

where  4>  is  the  prescribed  equivalence  ratio. 

2.  Compute  initial  values  of  P^'  where 

P,  =  KjK/’P4  %  * 
similarly  for  ?2',  P3*.  and  P5'. 

3.  Compute  initial  values  of  the  coefficients  of  the  8*  Eqs.  (11)  where 


=  0.5 P j/P484  -  0.254-5, 


similarly  for  82,  83,  and  85. 


5 


AEDC -TR-71-256 


4.  Compute  aj  and  a2  using  Eqs.  (12). 

5.  Compute  the  (9a/ SP^)  5^  terms  using  Eqs.  (13)  and  collect  all  the  coefficients 
of  64  and  56  together  to  form  a  2  by  2  matrix 


~°1  =  p^tO-5  Pj  +  0.5(l-rHQ)P3  +  ^-rH0)P4  +  2P5]S4 
4 

+  p-r[-0.2oP]  -  0.5rHOP2-0.25(l  — r„0)P3  +  P5  -  2rHOP6]56 


(14) 


~a2  =-p[o.5P,-o.5  p;+p;+pj*4 

4 

-*--i-[-0.25Pj  +  0.5  p'-r0.25P3  -  0.5  P5  +  Pg]S6 


(15) 


Substituting  for  56  from  Eq.  (6)  into  Eq.  (5)  gives  an  equation  for  64. 
Substitution  of  64  into  either  Eq.  (14)  or  (15)  gives  an  equation  for  56- 

6.  Compute  P4  =  p'  +  S4  and  P6  =  Pg  +  Sg 

7.  Repeat  processes  2  through  6  until  64  and  65  are  very  small  ( 1 0_s  is 
considered  sufficiently  small  for  most  calculations).  The  resulting  errors  in 
dependent  species  concentrations  are  always  less  than  this  (see  Eqs.  (11». 

8.  Compute  final  values  of  Pi,  P2,  P3,  and  P5  to  obtain  the  equilibrium 
composition. 

FUEL  RICH  OR  STOICHIOMETRIC  CASES 


P1  -  k,p5- 


p2  .  k2k;'p4p;* 


S1  =  0.5  —  ss 
5 

P„  P„ 

s2  =  Fs4--r«5 
•1  5 


p3  =  K3P4P£4 


P  P 

S3  =  /S4  "  °'5/S5 
4  5 


(16) 


Thus 


P6  =  K42P4P52 


=  2^54  -  2^S5 

4  5 


-al  =  ~'^~rHOP2  +  ^1-I'irn^P3  (2-iWP4  “  4rHOP6^6 


+  —  [O.o  P  ,  +  rFt0P  2  *  °-5d  “  rFI0^P3  +  2P5+  4rIIOP6^5 


(17) 


-  a2  -  -JP2  +  P3  +  P4  +  2Pg]  S4  +  ,[0.5  P  j  -  P2  -  0.5  P3  +  P5  -  2Pg]Sg 
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Solve  as  for  lean  case  with 

P5  =  P  -  P4  </>  >  1 

P'  =  10*2  <f>  =  1 

1.2  C/H/O/N/A  <m  =  3) 

In  this  example,  a  system  of  equations  is  derived  which  can  be  used  to  calculate 
equilibrium  composition  for  hydrogen/air  or  hydrocarbon/air  combustion  with  pressure 
specified  and  either  temperature,,  enthalpy,  or  entropy  specified.  This  set  has  been 
programmed  in  FORTRAN  IV  for  general  use,  and  a  copy  of  the  listing  is  attached  as 
Appendix  II. 

The  ‘following  n  dissociation  equilibrium  and  conservation  equations  completely 
describe  the  relationship  between  the  n  unknown  partial  pressures. 


i 

Fuel  Lean 

Fuel  Rich 

1 

P1 

y  y 

V  D'ip'4 

=  R.jr8  r? 

P1 

-  K,P^ 

2 

P2 

y 

_  if  p'i 

P2 

-  K2P8P6"‘ 

3 

P3 

y  y 

v  p'ip'« 

=  l'‘3r8  r7 

P3 

-  w;1 

4 

P4 

y  y 

V  Ti'l  Ti'l 

=  ^4^10^7 

P4 

y  y  -1 
-  K.P8,PWP6 

5  P5  =  K5P9P^  P5  =  KSP9P6P-B1 

6  p6  =  k6p8p;1/2  p7  =  k6p^p;2 


(18) 


Species 

1  =  123  4  5  67  8  9  10  Ar 

H  0  OH  H20  NO  CO  H2  02  C02  N2  A 

Conservation  Equations 
Lean 


P9=  rCN(2P10  +  P4)/(l  +  K5P7J/3) 
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Rich 


.P9  =  rCN(2P10  +  P4)/(l  +  K5P6/P8) 

PAR  =  rAN^2P10  +  P4^ 

P1  +  P3  +  2P6  +  2P8  +  rHN^2P10  +  P4^  =  ® 


P2  +  P3  +  P4  +  P5  +  2P7  +  P8  +  2P9  ”  rON^2P10  +  P4^  =  ® 
10 


(19) 


Epi  -<p-pa»>  - » 

1 

A  more  sophisticated  method  of  bookkeeping  is  used  in  order  to  simplify  the  program 
logic  in  this  example.  The  matrix  is  set  up  as  follows 


£  a-S<  =  a-  j  =  l»i 

i  ii  j  j  9 

i=l 

The  solutions  are  obtained  using  Cramer's  Rule 


(20) 


*1- 


-“1 

a12 

0 13 

~a2 

a22 

“23 

"a3 

a32 

a33 

(21) 


<aji) 


and  similarly  for  82  and  83 ,  where  J(aji)  is  the  determinant  of  the  coefficient  matrix. 


J(aij>  = 


11 

a  12 

a13 

21 

a22 

a23 

+  0 

(22) 

31 

a32 

a33 

For  fuel  lean  mixtures,  81  =  87 .  82  =  5s ,  and  83  =  81 0 !  for  fuel  rich  or  stoichiometric 
mixtures,  Sj  =  8g  with  82  —  8s  and  83  =810.  The  coefficients  in  the  matrix  equation 
are  derived  as  in  the  example  given  in  section  2.1.1. 


Oj  =  pi  +  P3  +  2Pg  +  2Pg  -  rHN(2p10  +  P4) 

a2  =  p2  +  p3  +  p4  +  p5  +  2P?  +  Pg  +  2Pg  -  rON(2P10  +  P4) 

10  t 

*3  =  Z  Pi  "  (P-PAR> 
i=l 


(23) 
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an  =  d^.  +  d3ti  +  2d6(i  +  2d8(i  -  rHN(2d10  -t  +  d4  j) 
a2i  =  d2,i  +  d3,i  *  d4,i  +  d5,i  *  2d7,  i  +  d8,  i  +  d9,  i 
10 

a3,  i  =  £  dK,  i  ‘r  dAR,  i 
K=  1 


rON^2d10,i  + 


d4  i> 

4’ '  (23) 

(Continued) 


FUEL  LEAN  CASE 


du  = 

-0.25  Pj/P7 

dl  ,2  =  Pl'/P8 

O 

II 

CO 

'O 

d2, 1  = 

0.5  P/P-' 

d2,2  =  ® 

d2,3  =  0 

d3 , 1  = 

0.25  P3/P7 

d3,2  =  P3^P8 

d3,3  =  0 

II 

no 

0.5  P 4/ P 7^ 

d4,2  =  0 

d4,3  =  °-5FVP10 

d9, 1  = 

rCN^d4, 1  *  °-5t2P 

;0+pX/p;  S/d 

+  K-/P7  )]/(1  +  Ks/P. 

d9,2  = 

0 

d9 ,3  =  rCN^2  +  d4,3 

;)/(l  +  K5/P7  S 

d5,l  = 

/  / 

p-  Pc 

d9  1—  -  0-5  — 

P  P 

9  r7 

dS,2  =  0 

d  d 

d5,3  =  d9,3  — 

P9 

d6.,  - 

-0.5  Pg/P, 

d6,2  =  P6/P8 

d6 ,3  =  0 

d7.1  = 

1 

d7,2  =  0 

d7 .3  -  0 

d8, 1  = 

0 

d8,2  =  ^ 

d8,3  =  0 

d10,l  = 

0 

d  10,2  =  0 

d  10 ,3  =  1 

dAR,  1  “ 

rAN  d4, 1 

d  A  R ,  2  =  0 

dAR,3  =  rAN(2+d4,3' 

FUEL  RICH  CASE 

dU 

=  0-5  P y'K 

d  1 , 2  =  0 

d 1 ,3  =  0 

d  2 , 1 

=  "P2/P6 

d2,2  =  P2,/P8 

d2,3  =  0 

d3, 1 

=  di.|P>; 

d3,2  =  P3/P8 

d3 ,3  =  0 

d4, 1 

/  * 

=  ~  'VP6 

d4, 2  =  P4/P8 

d4,3  =  °-5P4/ 

d7, 1 

=  -2P>6 

d7 ,2  =  2P7/P8 

d7,3  =  0 

y  (24) 


(25) 
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d9pl  =  ^rCNd4,l  “  K5P8P9^1  +  K5P6/,P8^ 

d9p2  =  ^rCNd4,2  *  K5P6P9,/P8  ^  +  K5P6'/P8^ 

*9,3  =  rCN^2  +  ^.a*1  +  K5P6/P8^ 

dS, 1  =  (d9,lK5/P8+P5)/P6 

d5,2  “  (d9,2^5P6  “  P5^P8 

d5,3  =  (dg^KsP^/Pa 


(25) 

(Continued) 


d6,l  = 

1 

d6,2  =  0 

d6,3  " 

0 

11 

•-H 

CO 

"T3 

0 

d8,2  =  1 

d8,3  = 

0 

d10,l  = 

0 

d10,2  =  0 

d10,3  = 

1 

dAR,  1  = 

rANd4, 1 

dAB,2  =  rANd4,2 

dAR,3  = 

rAN(2  +  d4>3) 

METHOD  OF  SOLUTION 

1.  Assume  values  for  Pg',  Pi0,  and  P7'  or  Pg';  evaluate  rAN,  rcN , rn n >  and 
Ton  using  the  stoichiometric  equation 


xC  +  yH  +  zO  +  aN  +  bA 

xC02  +  JL  H20  +  (-|-x-  -I)02  +  2  N2  +  bA  (Lean)  (26) 

or 

-j  CO  +  CO 2  +  (z-1.5x)H20  +  (0.5y  -  z  —  1.5x)  H2  +  -j  N2  +  bA  (Rich) 

Rich  Fuel  Case 

P6'  =  P(-|  -  z  +  1.5x)/[x  +  0.5(y+a)J 

Pg  =  P(z  -  1.5  x)/[x  +  0.5(y  +  a)]  (27) 

Pj0  =  P  a/[x  +  0.5(y  +  a)] 

Lean  Fuel  Case 

P?  =  P  [z  -  (2x  +  0.5y)]/(0.5  y  +  z  +  a) 

p  g  =  Py/(0.5y  +  z+ a)  (28) 

Pjq  =  Pa/(0.5y  +  z  +  a) 
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Atomic  Ratios 

rAN  =  b/a  rCN  =  c/a  rHN  =  >‘/a  r0N  =  z/a 

2.  Compute  P;'  using  the  dissociation  equilibrium  (Eqs.  (18)) 

3.  Compute  djj,  ajj,  and  dj  using  Eqs.  (23)  and  (24) 

4.  Solve  for  5j  as  in  Eq.  (21). 

5.  Correct  Pg Pio',  and  either  Pg’  or  P7'  using  the  6j  (Pj  =  Pj'  +  6j). 

6.  Compute  new  Pj',  i  =  1,  5,  P7',  and  Pg',  using  the  dissociation  equilibria. 

7.  Repeat  cycle  3  through  6  until  the  6j  are  less  than  a  specified  magnitude 
(10- 5  is  usually  adequate).  If  there  is  no  carbon  in  the  fuel,  then  x  should 
be  set  equal  to  zero.  If  there  is  no  carbon  in  the  system,  tcn  is  zero  and 
the  "lean"  ds  i  and  dg  3  calculations  should  be  bypassed. 

SOOT  FORMATION 


If  there  is  insufficient  oxygen  to  bum  all  the  carbon,  then  soot  will  be  produced. 
This  case  is  handled  by  adding  the  carbon  formation  equation 


Cs  +  C02  ^  2 CO  Ks  =  Pg/P \  <29> 

The  carbon  vapor  (saturation)  pressure  can  be  obtained  using  Gibb's  phase  rule, 
assuming  coexistence  of  condensed  phase  and  vapor  phase 


Pv  = 


*pKF 


„  ° 
C,R 


-  FCiW)/HTI 


(30) 


where  Fc°  is  the  standard  free  energy  of  carbon  and  subscripts  s  and  v  denote  solid  (soot) 
and  vapor  phases,  respectively.  Actually,  Pv  *  0  for  most  applications;  thus,  only  Ps  (soot 
concentration  -  not  a  pressure)  needs  to  be  considered  and  it  is  evaluated  using  the  carbon 
conservation  equation  (thus  P9  and  P5  are  modified  to  include  Ps).  If  Pv  is  not 
approximately  zero  (T  >  4000°C)  then  it  is  also  added  to  the  carbon  conservation  equation 
and,  in  addition,  is  added  to  the  sum  of  the  partial  pressures  (Ps  is  not). 


2.2  CALCULATIONS  WHEN  TEMPERATURE  IS  UNKNOWN 


2.2.1  Given  Static  Enthalpy  or  Total  Enthalpy  and  Velocity 

The  static  (sensible  plus  chemical)  enthalpy  (H)  is  calculated  from  the  current  values 
of  partial  pressure  and  the  initial  or  current  estimate  of  static  pressures  and  the  initial 
or  current  estimate  of  static  temperature  (T).  If  the  given  value  is  denoted  by  H*,  then 
the  correction  equation  for  temperature  is  obtained  from  the  enthalpy  equation 
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n 


i=l 


(31) 


where 

n 

hi  =  ■/’Cp.dT  T  Ahform  i  andPw  =  £  PiWi  (32) 

i=l 


P(HdW  +  WdH) 


(33) 


For  temperatures  below  about  4500“R,  the  frozen  specific  heat  is  essentially  equal 
to  the  equilibrium  specific  heat,  i.e., 


r  ^L 

**  dp. 
i=l  1 


dP.  -  PCTdW  «  dT 
1  d  1 


(34) 


thus  Eq.  (33)  is  approximated  quite  well  (for  calculation  purposes)  by 


or 


H*  -  H  =  AH  «  CpAT 
AT  -  AH/Cp 


(35) 


(36) 


The  corrected  temperature  T  =  AT  +  T. 

Using  Eq.  (35)  results  in  a  considerable  saving  in  the  number  of  calculations  which 
must  be  made  and  does  not  require  an  increase  in  the  order  of  the  matrix  which  must 
be  solved.  The  converged  solution  is  the  same  whether  Eq.  (34)  or  (35)  is  used  for  the 
correction  equation. 

2.2.2  Given  Entropy 


The  entropy  is  used  as  a  control  variable  in  isentropic  expansions  in  which  it  is  desired 
to  obtain  the  temperature  and  composition  as  the  pressure  is  reduced  during  an  expansion. 
The  correction  equation  used  is 


* 

S  -  S  = 


(37) 


or 


AT  =  T/C  AS 


(38) 
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The  corrected  temperature  T  =  T  +  AT. 

Above  T  =  4500°R  the  difference  between  the  frozen  and  equilibrium  values  of  specific 
heat  becomes  significant  and  difficulties  are  experienced  in  obtaining  convergence  of  the 
iteration  equations.  Some  additional  logic  controls  were  incorporated  in  the  computer 
program  to  prevent  stable  oscillations.  This  technique  has  proved  to  be  satisfactory  for 
temperatures  below  5500°R. 


SECTION  III 
COMPUTER  PROGRAM 

A  computer  program  was  written  for  the  Raytheon  520  digital  computer  to  solve 
for  the  equilibrium  composition  and  thermodynamic  properties  of  a  chemical  system 
containing  carbon,  hydrogen,  oxygen,  nitrogen,  and  argon,  specifically  hydrocarbon  air 
combustion  systems.  However,  because  of  the  versatility  of  the  input,  other  types  of 
problems  may  be  solved.  These  include  most  fuel  oxidizer  combinations  which  contain 
at  least  hydrogen,  oxygen,  and  nitrogen  in  the  fuel  and  oxidizer.  For  example, 

1.  Hydrogen  or  hydrocarbons  with  oxygen-enriched  or  vitiated  air 

2.  Hydrazine  or  methyl  hydrazines  with  N2O4  or  nitric  acid 

There  have  been  some  difficulties  obtaining  a  solution  to  some  types  of  problems.  These 
are  discussed  near  the  end  of  this  section.  The  program  listing  is  presented  in  Appendix 
II  and  the  results  of  some  calculations  are  given  in  Appendix  Ill. 

The  computation  was  divided  between  two  subroutines  which  are  called  by  a  main 
program  whose  function  is  to  input  and  output  data.  The  main  program  will  not  be 
discussed  in  detail  since  it  is  assumed  that  any  other  user  would  choose  to  rewrite  this 
section  for  his  own  specific  purpose.  The  equilibrium  calculations  are  made  in  subroutine 
PROP  which  calls  on  subroutine  THERM  for  the  required  thermodynamic  properties. 

There  are  four  optional  calculations  which  can  be  made  using  PROP.  These  are 
controlled  by  the  integer  variable  IT  in  the  argument  list,  and  are  discussed  below. 

IT  =  0.  Given  the  moles  of  the  elements  carbon  (CAR),  hydrogen,  (HYD),  oxygen 
(OX),  nitrogen  (AN2)  and  Argon  (AR),  the  pressure  (PR)  in  atmospheres,  and 
the  enthalpy  (H)  in  Btu/lb,  plus  an  initial  estimate  for  the  temperature  (TT) 
in  °R,  the  subroutine  returns  the  equilibrium  temperature  (TT)  in  °R,  the  entropy 
(S)  in  Btu/lb°R,  the  specific  heat  (CP)  in  Btu/lb°R,  the  molecular  weight  (WTM), 
an  array  of  the  partial  pressures  of  the  1 0  species  (P)  in  atmospheres,  the  partial 
pressure  of  argon  (PAR)  in  atmospheres,  and  the  number  of  iterations  (NT) 
required  to  obtain  the  solution. 

IT  =  1.  Given  the  moles  of  the  elements,  the  temperature  (TT),  and  the  pressure 
(PR),  the  subroutine  returns  the  enthalpy  (H),  partial  pressures,  and  the  other 
properties  returned  when  IT  =  0. 
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IT  =  2.  Given  the  moles  of  the  elements,  the  entropy  (S),  and  the  pressure 
(PR),  the  subroutine  returns  the  composition,  temperature,  and  other 
thermodynamic  properties  returned  when  IT  =  0. 

IT  =  3.  Given  the  composition  (array  of  P's  and  PAR),  the  entropy,  and  the 
pressure,  the  subroutine  returns  the  thermodynamic  properties  for  a  frozen 
composition.  In  this  case,  if  a  series  of  entries  to  the  subroutine  are  made  at 
different  pressures  (for  example,  a  frozen  isen  tropic  expansion  to  a  series  of 
known  pressures)  the  previous  composition  is  stored  so  only  the  initial 
compositions  must  be  entered. 

Normally  for  this  case  and  the  previous  one  (IT  =  2)  the  input  (S)  and  compositions 
are  obtained  by  running  a  type  0  or  1  case  immediately  before  the  type  2  or  3  series. 
For  example,  in  an  isentropic  nozzle  expansion  the  combustion  chamber  conditions  are 
first  calculated  from  the  fuel-to-air  ratio,  enthalpy,  and  pressure,  then,  using  the  entropy 
from  this  calculation,  constant  entropy  cases  (IT  =  2  for  equilibrium  or  IT  =  3  for  frozen 
expansion)  may  be  calculated. 

In  summary  the  argument  list  of  PROP  is  as  follows:  Inputs  which  are  unchanged 
are 


CAR  Relative  moles  of  carbon 

HYD  Relative  moles  of  hydrogen  atom 

OX  Relative  moles  of  oxygen  atom 

AN2  Relative  moles  of  nitrogen  atom 

AR  Relative  moles  of  argon 

PR  System  pressure  in  atmosphere 

IT  Type  of  case 

Variables  which  are  both  input  and  output  are 

H  Enthalpy,  Btu/lb.  For  IT  =  0,  H  is  entered  as  the  desired  mixture 

enthalpy  and  returned  as  the  enthalpy  of  the  mixture  found  by  the 
program.  For  all  other  cases,  it  is  an  output  only  and  is  the  mixture 
enthalpy. 

TT  Temperature,  °R.  For  IT  =  1,  it  is  the  desired  mixture  temperature 
and  is  not  changed.  For  all  other  cases,  it  is  entered  as  an  initial 
estimate  of  mixture  temperature  and  is  returned  as  the  converged 
mixture  temperature. 
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S  Entropy,  Btu/lb°R.  For  IT  =  4,S  is  entered  as  the  desired  mixture 

entropy  and  is  returned  as  the  calculated  mixture  entropy.  For  IT 
=  1  and  0,  it  is  an  output  only  and  is  the  mixture  entropy.  • 

P(10)  An  array  of  species  partial  pressures,  atm.  For  IT  =  3,  it  is  an  input 
using  the  values  from  the  previous  calculation  and  is  returned  ratioed 
by  the  pressure  in  the  previous  calculation  to  that  in  the  present  one. 

For  all  other  cases,  it  is  an  output  only. 

PAR  Argon  partial  pressure,  atm.  It  is  used  the  same  way  as  the  other  partial 
pressures. 

Variables  which  are  output  only  are 

CP  Specific  heat,  Btu/lb°R 

WTM  Molecular  weight 

NT  Number  of  iterations  required  for  convergence 

A  brief  description  of  the  subroutine  PROP  follows.  In  this  subroutine  the  equations 
are  solved  as  outlined  in  Section  2.1.2. 

1.  A  decision  is  made  as  to  whether  the  mixture  is  rich  or  lean  and  Initial 
estimates  are  made  of  the  partial  pressures  of  the  control  species  (H2 ,  H2  O, 
and  N2  for  a  rich  case  and  O2,  H2O  ,  and  N2  for  a  lean  case). 

2.  Inside  the  main  iteration  loop,  the  partial  pressures  (P(I))  of  the  remaining 
species  are  calculated  from  equilibrium  constants  calculated  from  free 
energies  obtained  from  THERM. 

3.  The  partial  derivatives  (D1(J)  to  DAR(J))  of  the  species  with  respect  to 
the  control  species  are  calculated. 

4.  The  coefficients  A1(I)  of  the  error  terms  DELI,  DEL2,  and  DEL3  and 
the  constants  SP1  and  SP3  in  the  expanded  conservation  equations  are 
determined.  The  matrix  is  solved  for  the  three  error  terms  which  are  then 
added  to  their  respective  estimates  of  partial  pressure  for  the  control  species. 

5.  If  T  is  known  (IT  =  1),  these  new  pressures  are  used  as  estimates  and 
the  process  is  repeated  until  convergence  occurs,  i.e.,  the  absolute  value 
of  DELI,  DEL2,  and  DEL3  is  less  than  DELL  (lO6  in  this  version). 

6.  For  the  other  cases,  an  enthalpy  or  entropy  is  calculated  using  the  new 
partial  pressures  and  the  estimated  temperature  by  calling  THERM.  The 
temperature  estimate  is  then  corrected  and  the  entire  process  repeated  with 
the  new  T  and  control  specie  pressures  until  convergence  occurs,  i.e.,  the 
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absolute  value  of  the  DEL's  less  than  DELL  and  the  error  in  H  dr  S  less 
th^n  the  frozen  value  of  specific  heat  or  3S/9T,  respectively.  This  control 
gives  a  maximum  error  in  temperature  of  1°K.  Subroutine  THERM 
computes  the  thermodynamic  properties  of  mixtures  of  the  eleven  species 
using  the  curve  fit  coefficients  from  Ref.  6.  It  is  called  through  an  argument 
list  consisting  of  the  following  terns: 

Inputs  not  changed  by  the  subroutine  are 

P  Array  of  partial  pressures  of  the  species 

T  Temperature,  °K 

PR  System  pressure,  atm 

PAR  Argon  partial  pressure,  atm 

KT  Integer  variable 

Outputs  are 

KT  =  0 

H  =  Mixture  enthalpy,  cal/gm 
CP  =  Mixture  specific  heat,  cal/gm°K 
F  =  Array  of  H/RT  for  10  species 
WTM  =  System  molecular  weight 
KT  =  1 

F  =  Array  of  F/RT  for  10  species 
H  and  CP  and  WTM  are  meaningless 
NOTE:  The  only  input  necessary  with  this  option  is  T. 

KT  =  2 

H  =  Mixture  entropy,  cal/gm°K 
CP  =  Mixture  S/T,  cal/gm°K2 
F  =  Array  of  S/R  for  the  10  species 
WTM  =  System  molecular  weight 
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A  sample  calculation  typical  of  that  done  in  a  main  program  is  outlined  in  detail 
in  Appendix  III. 

In  solving  a  combustion  problem  at  known  enthalpy  or  entropy  for  composition  and 
thermodynamic  properties,  there  are  four  unknowns.  These  are  the  three  control  species 
compositions  and  the  temperature.  The  error  terms  in  the  species  equations  neglect  the 
partial  derivatives  of  composition  with  respect  to  temperature,  and  the  error  term  in  the 
energy  equation  neglects  the  partial  derivatives  of  enthalpy  with  respect  to  composition. 
This  is  a  great  simplification  since  the  neglected  partials  require  time-consuming 
computation  and  also  add  an  additional  term  to  the  matrix  of  error  terms. 

In  some  test  problems,  the  solution  instead  of  converging  has  oscillated  around  the 
correct  answer.  These  have  been  cases  where  the  temperature  was  high  and  the  uncoupling 
of  the  species  and  energy  equations  was  potentially  serious.  No  cases  of  nonconvergence 
were  seen  when  the  final  temperature  was  below  5000°R  which  should  include  almost 
all  problems  of  interest  involving  air  as  an  oxidizer. 

Several  controls  have  been  added  to  the  program  to  increase  the  numerical  stability 
for  type  0,  2,  and  3  problems. 

1.  For  equivalence  ratios  near  stoichiometric,  the  initial  estimate  of  species 
concentration  for  H2  in  rich  mixtures  and  stoichiometric  or  Oz  in  lean 
mixtures  is  set  equal  to  0.01.  If  this  is  not  done,  the  initial  estimate  of 
some  of  the  free  radicals  will  be  abnormally  high  and  the  calculation  will 
not  converge. 

2.  When  the  error  in  temperature  changes  sign,  two  iterations  of  composition 
are  made  at  the  same  temperature.  This  control  has  prevented  the  solution 
from  oscillating  in  most  cases. 

3.  The  initial  estimate  of  temperature  has  been  left  as  an  input.  A  good 
estimate  will  increase  the  speed  of  convergence.  However,  it  is  usually  safer 
to  estimate  low  rather  than  high  since  the  equilibrium  constants  are  set 
up  to  be  large  at  high  temperature  and  thus  can  greatly  overpredict  the 
dependent  species  if  the  temperature  estimate  is  too  high. 

4.  The  change  in  temperature  per  iteration  step  is  reduced  below  that 
calculated  by  the  enthalpy  equation. 

SECTION  IV 

RESULTS  OF  CALCULATIONS 

Some  results  of  calculations  obtained  using  this  program  are  given  in  Appendix  III. 
Table  I  is  a  computer  printout  listing  the  results  from  a  number  of  cases  where  equivalence 
ratio,  enthalpy,  and  pressure  were  specified  for  a  CnH2n  hydrocarbon  burning  in  air.  Table 
II  lists  the  results  of  cases  where  equivalence  ratio  temperature  and  pressure  were  specified. 
Table  III  lists  the  results  of  a  frozen  entropic  expansion,  and  Table  IV  lists  the  results 
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of  an  equilibrium  isentropic  expansion  from  10  atm  to  a  given  pressure  for  an  equilibrium 
gas  with  a  given  equivalence  ratio  and  temperature.  These  results  have  been  compared 
with  those  of  Ref.  7  and  computed  results  using  the  program  of  Ref.  1.  For  the  same 
input  conditions  the  temperatures  agree  to  within  1°K,  the  compositions  of  the  major 
species  to  less  than  1  percent,  and  the  minor  species  to  a  few  percent. 

The  number  of .  iterations  required  shows  that  problems  where  composition  or 
temperature  only  is  a  variable  (frozen  isentropic  expansions  or  known  temperature) 
converge  more  quickly  than  problems  where  both  temperature  and  composition  are  varied. 
It  also  shows  that  solutions  are  obtained  more  quickly  at  lower  temperatures  where 
dissociation  is  small. 

The  times  required  for  the  calculations  shown  in  Table  I  are  for  the  Raytheon  520 
computer  and  would  be  expected  to  decrease  if  a  faster  machine  were  used  (approximately 
a  factor  of  five  for  the  IBM  360). 
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APPENDIX  I 

THERMODYNAMICS  OF  CHEMICAL  EQUILIBRIUM 
MIXTURES  OF  PERFECT  GASES 

A  chemical  reaction  is  formulated  generally  as  follows: 


(p/ij  =  (P.pi 


The  free  energy  per  mole  of  a  mixture  is  given  by 

n  P.  n  p. 


f  -I  T  F1  -  E  T^T.i  ♦  RTJ"p|) 


i=  1 


i=l 


The  logarithm  of  the  equilibrium  constant  for  the  reaction  is  given  by 


f"Ki  -  kJjlPjj)  1J  ,J  =  £  -  ^j)«"  Pjj 

i=l 

The  change  in  free  energy  in  a  chemical  reaction  is  given  by 
iFi  .£ 


i=] 

n 


//  / 


-ES-'>p(fT,.  +  nTt"[y 


i=i 


//  .  o 


=  RTfnKj  (vVi  -  *.,)  rT  i 

i=i 


Chemical  equilibrium  is  readily  defined  as  follows: 

AFj  =  0 
K  =  exp 

All  enthalpies  for  the  reference  states  of  the  elements  are  0  at  77°F  (298.15°K). 


Hfuel  is  the  molar  enthalpy  of  the  fuel  at  the  inlet  temperature.  For  a  compound 
such  as  the  fuel,  or  the  CO2  in  the  air,  or  an  element  not  in  the  standard  state,  the 
molar  enthalpy  is: 

T,  °R 


I 


53(5.67 


CpdT 


H 


f  o 

536.67  R 
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for  example,  air  at  540°R 

H  =  0.780881  x  23.22  =  18.132 

+0.209495  x  23  =  4.902 

+0.909324  x  16.56  =  0.154 


+0.00030  x  -169263.7  =  -50.779 

-27.590  Btu/lb  mole 

Heats  of  formation  are  not  available  for  most  mixed  hydrocarbon  fuels  such  as  JP4. 
The  energy  content  of  this  type  of  fuel  is  more  commonly  given  as  heat  of  combustion. 
The  relation  between  these  two  properties  may  be  derived  as  follows: 

^reactants  (T  ref)  =  ^products  (X  products) 

"  Qcomb.(T  ref)  Hpro<|ucti  (T  ref) 

The  reaction  for  a  CH„  fuel  is 

CHn  +  (1+  f)02  -  C02  +  1  H20 

for  a  basis  of  1  mole  of  fuel  with  a  heat  of  combustion  of.  Qcomb.  x  WTMfuei 


H fu el  (T  fuel)  +  Cpfuel  x  (Tref-Tf«eP  +  (1_f)CP0  (T~77)  =  Qcomb.(Tre() 
x  (12.011  +  1.008n)  +  Hco^(t  ref)  +  -  H{|20(T  ref) 


For  example,  a  CH2  hydrocarbon  with  a  lower  heating  value  of  18,000  Btu/lb  at  70°F 
and  a  specific  heat  of  0.53  Btu  would  have  an  enthalpy  at  20°F  calculated  as  follows: 


H(uel  (T=  20)  =  -°-5  x  14.027(70-20)  -  (l  +  x  7.02(70  -77) 
+  18000  x  14.027  -  169.335  -  104092 
=  -21410  Btu/lb  mole 


If  this  fuel  is  mixed  with  540°R  air  at  an  equivalence  ratio  of  0.5,  the  enthalpy  entered 
into  the  program  is: 


-27.59  +  0.5  x  -  21410 
28.9666  +  0.5  x  1.95905 

=  -358.4  Btu/lb 
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APPENDIX  II 

COMPUTER  PROGRAM  LISTINGS 
NOMENCLATURE  FOR  MAIN  PROGRAMS  INPUTS 

IT  =  0 

NPHI  =  Number  of  equivalence  ratios 
NP  =  Number  of  pressures 
NH  =  Number  of  enthalpies 
PHI  =  Equivalence  ratio 
P  =  Pressure,  atm 
H  =  Enthalpy,  Btu/lb 

IT  =  1 

same  as  IT  =  0,  except 
NH  =  Number  of  temperatures 
T.  =  Temperature,  °F 

IT  =  2  or  3 

same  as  IT  =  1,  except 
T  =  Stagnation  temperature,  °F 
ITT  =  IT 
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MAIN  PROGRAM  FOR  KNOWN  ENTHALPY,  IT  =  0 


DIMENSION  XXU4.50) 

_ 0 1  HEX'S  I  ON  PHI  ( 1 0 ) . P( 10 ) <  H (  lO).X(lO) _ 

READ  ( 29# 97 )  NPHI.NP.NH 

_ RE  A.Q.J2.9.t  96J_(  _ 

READ  (29,96)  (P< J).J«1.NP) 

_ R  E  AJLJ.2  £,.*6  1  ._NH ) _ _ 

99  FORMAT  <  2F12.4,F10.1,F12.2.2F12.5,F11.3,2I6.F9.3) 

_ 9.8..F0R»AT...U,Ml2t4i _ _ 


97  FORMAT  (315) 

96  FORMAT  (PElOol _ 

95  FORMAT  (8X, 2HE0, 10X, 1HP.9X. 1HT , 11X , 1HH , 11X , 1HS , 12X , 2HCP, 7X , 3HM0L , 

*5X,2MN0,_4Xj_4HT_vPE,2Xi4HTlM£.) _ 

94  FORMAT  (8X.5HRAT!0.6X,3HATM.6X,5HDEG  R ,  7X, 6HBTU/LB. 5X , 0HBTU/L8  OR 

_  *•  4X, 9HBTU/LB  1R . 4X , 2HWT , 5X. 4H I TE R , 3X. 4HC ASE, 3X , 3HSEC/ ) _ 

93  FORMAT  (5X.2HEQ  , 6X. IMP# 7X» 1HT, 12X » 25HSPEC IES  MOL  FRACTIONS) 

_  92.  FORMAT  (4X,5HRATI0.3.X,3iiAtM,4X.J.HD_EG  R.7X.1HH,  11X .  1H0 , 11X , 2H0H. 


•10X.2HN0. 10X,2HC0/> 

91  Format  (4X.5HRATI0,3X,3HATH,4X . SHDEG  R.7X, 2HH2.10X,2H02,9X _ 

•.3HH20.9X, 

»3hCQ2,10X,2HN2.10X,1HA/) _ 

90  FORMAT  ( 80H1  SPECIFY  ENTHALPY  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  A 

_ •  QH2  HYDROCARBON_H1TH_AIRI _ _ 

89  FORMAT  (F9.3.F7.2.F8.1.6E12.4) 

86  format  (F9.3,F7.2.F8.1,5E12.4) _ 

Ca.  41959 


r.  =  l. 561762 _ 

Ea. 009324 

KK  a  0 _ 

WRITE  (21,90) 

WR I TE  ( 21,95)  _ 

WRITE  (21.94) 

IJ^O _ 

DO  1  Jal.NPHI 
A  a . 1396633»PHI (J)-»,0003 
3a .279326*PHI ( J) 

DO  2  JJal.NP _ 

DO  3  JJJal.NH 


1  =  4000^ _ 

CALL  PROP  (A.B.C.D.E.T.P(JJ),H(UJJ).S,CP,WT,X,PAR.TIME.NT,IT) 

IF  (NT.GE.50)  GO  TO  2 _ 

WRITE'  (21,99)  PHI(J),  P(JJ),T,H(JJJ),S. CP. HT.NT.IT, TIME 


_  iKaKK  +  t _ _ j _ 

RPal./P(JJI 

_  DO  4  K«4, 13 _ _ 

4  X X ( K  ,KK)=X(K“3)*RP 

XX  (1  ,  KK  )  aPH  I  (  J ) _  _  _ 

XX ( 2  .KK)aP(JJ) 

__  JLXJ  3  .  KK)a  T _ 

XX ( 14 «KK)=PAR*RP 

._3_C0N_T INUE _ 

2  CONTINUE 

1  CONTINUE _  _____  _ 

WRITE  (21.90) 

WRITE  ( 2l'.  93) 

_  WRITE  (21.92) _ _ . _ 

"write'  (21,88)  ( (XX(J.K) » J*I, 8) . K a  1 , KK ) 

WRITF  (21.90) _ 

WRITE  (21.93) 

WRITE  (21,91)  _  _  _  _ 

WRITE  (21,89)  <(XX(J»K)#Jsl.3),(XX(J,K>.J=  9 , 14 ) ,K»1» KK ) 

STOP  .  _ _ 

END 
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MAIN  PROGRAM  WHERE  TEMPERATURE  IS  KNOWN.  IT  =  1 


01  MENS I ON  XX(14,'50> 

DIMENSION  PHI  (1<1)  ,P(10>  ,T{  J  0>  .X<10> 

MEAD  ( ?9 .97 )  NPHI.NP.NH 

MEAD  (29,96)  (PriH J),Jal.NPHI  ) _ 

MEAD  (29.96)  (P( J) . J=1,NP) 

READ  (29.96)  (T(J)jJal  .  NH) 

99  FORMAL  '(  2F12 . 4 • F10 .1 »F12 .  2  j  2F12 , 5 , Fll . 3V2 1 67KT73T  " 

93  FORMAT  (1CE12.4) _  _ _  _ 

'97  FORMAT  (315)'  .  . . 

96  format  ( PElfl .  ) _ _ _ 

95  FORMAT  (8X.2WE0*nX.lHP,9X.lHT,llX.lHH,UX»lHS.12X.2HCP#7X.3HM0L. 
•5X.2HN0.4X.  4«T vpt , 2X .  4HT  IME ) 

94  Format  (*X,5HRAT  I0.'6x;'3lTrTM,6X'7'5h''D'E''G  R7?5G6H0TU/LB75X,9R9TU/L0'  0R“ 

..  4X  .  9HQTU/LB  OR , 4 X .2HHT, 5X . 4H I T ER , 3X , 4HC ASE . 3X . 3HSEC/ ) 

93  FORMAT  “(5X  .2HE0  .  6X .  1HP ,  Tx.  iHTTl2X,'25HSPEC IES“  MOL  '  TRACTIONS) . 

92  FORMAT  (AX.SHRATIO.SX.SHATM.AX.SHDEG  R ,  7  X , 1HH .  11X . 1H0. 11X . 2HQH. 

*10X. 2HNO.10X.2HCO/  ) 

91  FORMAT  (4X.5HRATI0j3Xj3HATMi4X.5HDEG  R. 7X . 2HH2 . 10 X. 2H02. 9X 
•  » 3HH20 , 9X  . 

•3HCO2.10X,2HN2,i0X.lHAA) 

9b  Format  («3hi  specify  temperTuFTEV  equivalence  ratio  and  pressure  fo" 

*R  A  CH2  HYDROCARBON  WITH  AIR> 

8$  FORMA  T  (  F9 . 3»  F7 . 2i  FSTT .  6El2 . 4  > 

88  FORMAT  (F9.3.F7.2.F8.1.5E12.4) 

C=. 41959 
Dal. 561762 
E=.  009 '3  2 '4"“ 

4K  =  0 

— wR,rE  <2i;w5 - 

HRITF  (21.95) 

wRTIT  (21,94)  "  : 

I  T  =  1 

no  i  j=i,nPhi 

A=.1396633*PMI ( J)*.0003 

- 9a.T79326.PHlU) - 

00  2  JJ=1,NP 
no  3  JJJal.NH 

CALL  PRO°  (A,B,C.D.E,T<JJJ),P(JJ).H,S,CP,WT,X,PAR,TIME.NT,IT> 

- IT- CKHTGETTO)- GD— TO- 2 - - - 

WRITE  (21.99)  PHKJIi  P  C  JJ>  .  T  ( JJJ ) .  H.  S.  CP .  WT .  NT.  I T  ,  T I  ME 

<K=KK+1 

RPal./P(jJ) 

.  DO  '4  "Ki'4.r3  .  -----  - 

4  XX (K  ,KK)=X(K-3)*RP 
XX(1  .  KK  )  ap'Hl  (  J) 

XX ( 2  .KK)aP(JJ) 

XX (3  .KX)ST(JJJ) 

XX(l«.KK)aPAR*RP 
3  CONTINUE 

2  CONTINUE _  _ _ _ 

1  CONTINUE 

WRITE  (21.90) _ _ 

WRITE  (21.93) 

WRITE'  (21.92) 

_ WRITE  (21.88)  ((XX(J,K)»Jsl.6),K  =  l,KK) _ _ _ 

WRITE  ( 2i. 90) 

_ WRJ.TE_(_21  .93) _ _ _ _ 

WRITE  (21.91) 

WRITE  (21,89)  ((XX(J.M),Jal,3).(XX(J,K).Ja  9,14),Kel,KK) _ 

STOP 

END  ....  _ _  _ 


25 


AEDC-TR-71-266 


MAIN  PROGRAM  FOR  ISENTROPIC  EXPANSIONS,  IT  -  2  OR  3 


DIMENSION  XX ( 14, 50 ) 

_ DIMENSION— PiUXlJl.)  ,P(10 )« T  (10 ).  X(lflJ _ , - = 

100  CONTINUE 

- SEAD-Uft  »-fl  J-X-NP-H1  «i(P  «NN«JXI - 

IF  (NPHI.EQ.'IJ  GO  TO  101 

- fiEAB-12SUQ6J_lDHiXilW^J(gJiI-) - : - 

READ  (29*96)  (P< J) # J«1,NP) 

- READ.  (.2-9*  961—  (-T-UH * J*1 .NH ) - 

99  FORMAT  (  2F12 .4,F1Q.1,F12.2.2F12.S.F11.3,2I6.F9.3> 

96  FORMAT  tlflElg .41 _ 

97  FORMAT  (413) 

_9t6.  FORKAJL  4  &E10.-Q.1 _ _ * _ 

93  FORMAT  <8X.2HE0.10X,1HP»9X»1HT.11X.1HH.11X.1HS.12X.2HCP#7X.3HNOL. 

_ ftSX*2mtO.«.*X..4RTYP.E.2XJ6HTlHE) _ _ _ 

94  FORMAT  ( 8X.5HRATI0, 6X , 3HATM. 6X.5HDEG  R, 7X. 6HBTU/L8.5X, 9HBTU/LB  OR 

_ JL. _ 4Y.9HRTU/I  H  OR . 4X. 2HWT. 3X* 4H 1  TER, 3X, 4HCASE* 3X. 3HSEC/ ) _ 

93  FORMAT  (5X*2HE0  . 6X, 1HP* 7X*1HT, 12X . 23HSPECIES  MOL  FRACTIONS) 

•10X  *  2MN0 .  1  OX .  2MCO/ ) 

•  * 3HH20.  9X* 

_ t3HCQ2*10Xi2HN?-*.10X.jlHAAJ _ 

90  FORMAT  (83H1  SPECIFY  TEMPERTURE  EQUIVALENCE  RATIO  AND  PRESSURE  FO 

_ *&_ A -CH2 -HYDROCARBON  Nl.TH  AIR) _ 

89  FORMAT  (F9.3.F7.2.F8.1.6E12.4) 

FORM  AT  _IF9. 3,F7.2^6d.,3El2 .4 ) _ 

85  FORMAT  (45H0  FROZEN  ISENTROPIC  EXPANSION  FROM  10  ATN.) 

■J6  FORMAT  (1HQ) _ 

87  FORMAT  (45H0E0ULIBRIUN  ISENTROPIC  EXPANSION  FROM  10  ATM.) 

_ C-..41939 _ 

D«l, 361762 

_ Ej>*.0.Ql9124 _ 

KK*0 

_ HRJTE-A21ji.9DJ _ 

IF  (ITT.E0.2)  WRITE  (21.87) 

- LE--UJT.SQ  WRITE  —121x851 _ 

WRITE  (21.93) 

_ HR1TE-121j_9A)- _ 

DO  1  J»1 , NPH I 

_ <L«J.119it633«PHl(.J)t...Q(m3 _ 

B».279326*PHI( J) 

_ D0l_2  JJJf  IjJW _ 

TX*T( JJJ) 

_ Usl _ 

DO  3  JJ-l.NP 

_ LF_  (JJ-2)  IQ. 12.12 _ : _ 

12  IT«ITT 

_ JLX _ t_IX _ >(P(JJ)/P(JJ-1)  )»»(1.987^(CP#HT) ) _ 

10  CALL  PROP  (A.B.C.D.E.TX  , P( JJ) .H.S.CP, WT.X.PAR.TIME.NT, IT) 

_ IF.  (NT. BE. 50  )  00— 10— 2 _ 

WRITE  (21.99)  PHI ( J) «  P(JJ).TX  ,H, S. CP. NT, NT, IT, TIME 

_ 

RP»1./P( JJ) 
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ISENTROPIC  EXPANSIONS,  IT  =  2  OR  3  (Concluded) 


DO  4  K  =  4 , 13 

_ : _ 4.  XX(K.  ,KK)3X(K-3)«RE _ 

XX(1  « KK ) sPHI ( J  } 

_ -XX(2  ,KK) =P( JJ) _ 

XX ( 3  *KK)=TX 

_  XX(14,KK)aPAR»RP.  ...  _ _ _ _ 

3  CONTINUE 

_ _  2  CONTINUE  . - -  - - - - 

1  CONTINUE 

_ _  .  .  -WRITE  (21.90)-  ..  _ _ 

IF  UTT.EQ.2)  WRITE  (21,87)  ■ 

_ IF  (I.TTiEQ.3).  WRITE -(21*851 - 

WRITE  (21,931 

_ HR  I  TE- -1-21, 92  ) _ _  - - - 

WRITE  (21,88)  ((XX(J.K).J«1.B),K"1.KK) 

.  _  WRITE  (21,90) _ _ 

IF  (ITT.E0.3)  WRITE  (21,85) 

_ IF  (ITT. EQ. 2)  WRIT£_X21jJ2J _ 

WRITE  (21,93)  ' 

_ _  WRITE  _(2l,.91>_  _ _ 

WRITE  (21,89)  <<XX(j,K), J»l,3), (XX(J,K), J»  9, 14) ,K»1,KK) 

_  ..  GO.  TO  .100  _ _  _ 1 _ 

101  CONTINUE 

_ _ STOP  _  ...  _ _ 

END 
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SUBROUTINE  PROP 


SUBROUTINE  PR0P<CAR,HYD#0X,AN2,AR.TT,PR,H,S,CP.HTm.P»PaR 

•TIME. NT. IT) 

0 

c 

fi 

SPECIES  ARE  1-H, 2-0, 3-0H,4-N0,5-C0,6-H2, 7-02. 8-H20.9-C02 

SUBROUTINE  FOB  DETERM  I N I  NO  EQUILIBRIUM  PART  I AL  PRESSURES 

.10-N2 

c 

c 

T  IS  TEMP.  IN  DEG  R  AND  CAN  BE  EITHER  INPUT  OR  OUTPUT 
■ _ P  IS  PRESS  IN  ATMOSPHERES  ANO  IS  AN  INPUT _ 

c 

c 

H  IS  THE  ENTHALPY  IN  BTU  PER  LBM  AND  IS  EITHER  AN  INPUT  OR  OUTPUT 

S  IS  The  ENTROPY  IN  RTU  PER  LfiM-DEG  R  ANn  EITHER  INPUT  OR  OUTPUT 

c 

c 

CP  IS  THE  SPECIFIC  HEAT  IN  BTU  PER  LBM-DEG  R  AND  IS  AN  OUTPUT 

HTM  IS  THE  MOLECULAR  HT  OF  THE  MIXTURE  AN  IS  AND  OUTPUT _ 

c 

c 

TIME  IS  THE  PROGRAM  EXECUTION  TIME  IN  SECONDS  AND  [S  AN 
NT  IS  THE  NO  *  OF  ITERATIONS  TO  FIND  PARTIAL  PRESS  AND  IS 

OUTPUT 

AN  OUTPUT 

c 

c 

IT  IS  A  CONTROL  CONSTANT  FOR  THE  SUBROUTINE 

I Ta 0  H  AND  P  ARE  KNOWN  AND  T.CP  %  HTM  ARE  FOUND 

c 

c 

I T»1  T  *  P  ARE  KNOWN  AND  ARRAY  OF  PARTIAL  PRESSURES  ARE 
I Ta?  S  %  P  ARE  KNOWN  AND  T.CP.  s  WTM  ARE  FOUND 

FOUND 

c 

c 

I Ta3  T  *  P  »  PARTIAL  PRESSURES  ARE  KNOWN  -  FROZEN  CASE 

1 T«1  HAS  REEN  MODIFIED  TO  ALSO  CA1  CUI  ATF  M.S.fiP.WTM 

DIMENSION  D1<3),D2<3),D4<3),D5(3),D6(3),D7(3),D8(3), 

»(3).D3(3I.E(6)  . 

D9(3) »D10 

DIMENSION  DAR(3) 

DIMENSION  F(10],Ai(3),A2<3),A3<3),P(10>  .  _ 

C 

CHANGE  UNITS  OF  INPUTS  TO  BE  COMPATABLE  WITH  PROGRAM 
_ T  a  (5./9.)«TT _ : 

HO*  15. /9. >#H/. 99931 
__  _  _ IF  (_LT_tGEj_2_).  H0*S/. 999331 _ 

C 

IC.TtO  _  .  _ 

DELTL«0. 

CALL  TIMER  <TT1> 

IF  (IT  -  3)  3,2,3 

3  CONTINUE 

DELL". 000001 
_F_U£L»2.*CAR*.5«HYD 

C 

RATIO  OF  ELEMENTS 
_ RAN  ■ -AR/AN2 _ 

RCN*  CAR/AN2 

_ _ RHN»HYD/AN2  .  . 

.  c 

RON«OX/AN2 

DEL1«0 . 

DEL2«0. 

DEL3«0. 

IF(OX.GT.FUEL)  GO  TO  1 

IF  (CAR.LE.OX)  GO  TO  20 

WRITEC21.97) 

RETURN 

20  XXOj  /A  CAR*  .5*  (HYD  +  AN2)  1 

C 

INITIAL  GUESS  OF  H2  H20  AND  N2  PARTIAL  PRESSURES 

PJ.6  I  *PR*XX*  C .  5»HYD?0X*1_.5#CAR )  - 

P(8)»PR*XX«(  +OX-1 ,5*CAR ) 

PI10)bPR*XX«AN2*.5 

L'  6 

_ an  Tn  71 _ 

1  CONTINUE 
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SUBROUTINE  PROP  (Continued) 


C  INITIAL  GUESS  OF  02«  H20 . N2  PARTIAL  PRESSURES 

_ XX  5/ <  (  •  5«HYJ)*0XaAN2)JL.-5*AR-1 _ 

P(7)*PR*XX*(0X-FUEL) 

_ ELI  8)  aHYD*XX*.PR _ 

P(10>»  AN2*XX*PR 

7  G0_  T0~  71 

C  RESET  "PARTIAL  PARTIAL" PR E S SURES’ ~FOR  FR 0_Z eITc A SE--1T«3 

_ 2_C0NI1NUE  _  _ _ _ 

PXaPR/PHOLO 

_ D.g.72_  _ 

72  P(J)»P(J)«PX 

_ PAR.PPAR.aELX _ _ 

C 

_C _ ITE.RATI0N_.L0.0P.  10  .F_L.ND_P.ART I AL.  .PRESSURES _ 

71  CONTINUE 

_ LE-iP  CU.«LEa.4UJ_P.(UjjJ)A _ 

DO  4  K«1 » 50 

_ K0UNT*K _ _ _ _ 

C  SKIP  PARTIAL  PRESSURE  CALCULATION  FOR  FROZEN  CASE 

_ JF_  (IT  •  _3.)_.8j>_7.j_8 _ _ 

8  CONTINUE 

£ _ CHECK  ALL-JLN.LLLAL-PA.R Il'AL  £RES.SUHE5__EflR_ NEJSAl  I VE...  VALUES 

DO  41  J«l,10 

_ 4i_J£  J  P.< J_)  ..LI ,0.,J  _P_t.J.  E -3.  ...  _  _ 

C  DON'T  RECALCULATE  FREE  ENERGY  RELATION  FOR  T  KNOWN--IT«l 

_ .1  F.._C1. .  EQ  ».6_) _ GO.__T.D_-4  2 . . . _ _ 

IF  (K.GE.2.AND.IT.EQ.1)  00  TO  10 

_ CAkL-THJRH  (  RiLtH  ._5 P  tE i PR  r P_A.R_t_WHL._U _ 

E( 1 >  =  2.7182818285«»(.5*F(8)-.25«F(7)-F(1>) 

_ E12L 2... 7482 8JL82 8  5  •  •  t  A5  - F (2 ) )_  __  _ 

it  3)  a  2.7182818285#«( .5  *F f 8 > ♦ .25»F 1 7 > -Ft  3 > ) 

_ E  t  4 )  =  _2,_71828_18285»»t . 5*F { 10 >* . 5*F t 7 ) -F 1 4  LL _ 

Et5)  a  2 . 718 2818 285 •• f F ( 9) - . 5*F (7 } -F 1 5 )) 

_ E  Lfel_  a  2.7 1828182 * t  F(8)-,5 >F  t7)-FI6  )J _ 

C  CALCULATE  PARTIAL  PRESSURES  OF  THE  ASSUMED  KNOWN  SPECIES 

10  Ptl) ■_ _S Q RTj! P.t.8J. /SQRT(PC7) )  )«E (JJ _ 

P 1 2 )  *  SORT  t  P ( 7 ) > • E <  2 ) 

_P(3)aS0RTJP<8>«S0RTtP(7n  )#E(3)_ _ _  _ _ 

P(4)aS0RT(P(ir)«Pt7)>*E<4) 

_ XX » 1 .  /  (.  1,.*E  1  5  )  /SORT  tPt  7 )_>_) _ 

P ( 9 ) ■  RCN«t2.*P(10 >*P(4) )*XX 

_ P(5 L»_  PmzSORltJP (7_)  >*EL5_) _ _ _ 

P 1 6 ) ■  Pt  8 ) /SORT  <Pt7})*E(6) 

_ P.AS_*.J.2.i_JL_P_U.QJ  _*_P.t_4JJL*  R  AN _ 

C 

_ GO  TO  43 _ _ 

42  CONTINUE 

_ I F_iK lGEj 2.AND, 1T.EQ.1I  GO  TO  46 _ 

CALL  THERM  <  P»TiH*CP.F.PR»PAR»WTM»1) 

_ E 11 )  _«  _EJtP.JLi.5_*F.t  6.1-FJ  1 )  J _ _ 

E 1 2 ) B  EXPtF(8)-Ft6)-Fl2>) 

_ EilLi£XP_t£JL8.L=£J  lirJLUi  > _ 

E  1 4 } ■  EXP(F(8)*.5*F(10)-F(6>-F<4) ) 
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SUBROUTINE  PROP  (Continued) 


E<5>*EXP<F<6>>»F<9>-F<5>-F(8>  > 

_ £I6Jjl  _EXb_(  2-.*E-f-8 )  -  2  ,Jt£l  6L-E  L7A1 _ _ 

46  P  ( 1 ) >E ( 1 ) *S0RT <P(6 ) ) 

_ PtS>l«E(g1»P  IB  JZEiiJ _ 

P(3)aE<3)*P<8)/P<l) 

_ Pi4)tiet41«P(fl>»SQBTfPM0n/Pt6) _ 

P<7>»E(6)*P(B)*P(8)/<P(6)«P(6>) 

_ 3UT.«.E.t.5.L»P_(.6  _ : _ 

XX«1./(1.*YY) 

_ EiAURCN*L2  •.•PHIL)  *£_L4Jll*XX _ 

P(5)*P(9)*YY 

_ EA&-1 _ LZj-  »  FjJjaJL-aJElALieflAM _ 

43  CONTINUE 

_ C _ DERI  VAT  LON -OF.  J?ELAII.0NSM1P_0F.  KROJN  Se£CJ£JLWTWRESPECl_IO._QUESS„. 

DO  11  I al«  3 

_ D-l-U-U-tL. _ 

D2< I )»0, 

_ D3( I ) BQ . _ 

D4(I)»0. 

. . Q5.i.l ) .1.0  • _ 

D6( I)«0* 

_ D7U)a.lL. _ 

DB(I>-0. 

_ D9(M»0. _ 

11  DIO ( I ) *0 ■ 

_ LT^LL^JEQ._61_(m_TJ-44 _ _ . _ 

Dl(l)»  -PU)/P(7>*.25 

_ 01  i  2 1  «__P_U )  f?  ( 8J  • .  5 _ _ 

D2(l)a  P(2)/P(7)*,5 

_ D3(l)» _ P(3)/P(7I».23 _ 

D3<2)»  P(3)/P(B)*,5 

_ _ D4  (1 )  a  P(4I/P(7)a.5 _ 

D4(3)aP(4)/P(10)*.5 

_ IT_JJJA5j.E0tJljOJ_J30  TO. .441 _ : _ 1 _ 

D9 tl ) ■  RCN*(D4(l)-<:(2,»P(10)*P(4),)aXX*,5«El5)/  P(7)**1.5)*XX 

_ _ 

D5(l)a  P(5)/P(9)*D9<1)-P<5>/P(7)a.5 

_ D5(3)» _ P(g)/P(9)*D9<3)  _  _ ! _ 

441  CONTINUE 

_ Mil)  A )  /Pi  7 )  *_,  5_ _ 1_ _ 

D6<2)=  P(6)/P(8) 

_ ; _ 07ii  utuji _ : _ 

DB ( 2 ) =1 , 0 

_ D 1 0_(  3.1*  1...0 _ 

DAR(1>«RAN*D4C1> 

_ AAfUMsP* _ 

DAR<S)sR*N*<D4(3)+2. ) 

_ £D_  JflL.  4  5 _ _ _ 

44  CONTINUE 

_ pi  il)  a  p  Ul/P  _ _ : _ : _ 

D2(l)a-p<2)/p<6) 

_ D2  (  2)  e  P(g)/P(B) _ ; _ 

D3(l)a-P(3)/P(1)*D1<1) 

_ D1(2I«  P ( 3 ) /P ( B  ) _ 

04<l)a-P(4)/P(ft) 
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SUBROUTINE  PROP  (Continued) 


D4(2)=  P(4)/P(B> 

JLCJ  J./.R.UJUj.J. 
D7U)*-2.*P(7)/P<6> 


_ D7L21=_P_C.7J-/P.<8  )*2., _ 

IF  (CAR.EO.O.;-)  00  TO  442 

D9U)  =  XX*(RCN*D4<1)-P(9>*E<5>/P<8>) 

_ _D9(_3)*RCN*XX*  (2..  *04  (3)  ll _ 

D5<1)*YY*D9<1)*P<5>/P<6> 

_ 05  <  2  >  =  Y  Y_*D  9(2)  -P_<  5J  ZP_<  BJ _ _ 


D5(3)=YY*D9(3> 

_ 442  .CONTINUE _ 

D6(l)al. 

_ D8(2)=l, - - 

D10(3)»l. 

_ 0AR(L)  =  RAN*Q-4LU _ _ 

DAR(2)»RAN»D4(2) 

_ DARmsHAK*X2..i0i.(JJJ _ 1 _ 

45  CONTINUE 

c  “derivation  of'maTrTx"  coeftTci ents 

_ D012  I  *1  »_3_ .  _ 

A1(I)>  D1(I)+D3(])*2.*D6(I )*2.*D8( I ) -RHN*(2 . *010 ( I  )*D4( 1 ) ) 


_ A2CJJ-? — 02<L1*I 

**DlO<n*D4(I)) 

_ _ A3  <_I> «  .  01  <  I )  *02  U  )  *D.3L<  U  *D.4.LL>  *D5  (U  *D$1L)*P1.(  I  L*0JlJLI)*D9a)*DU9  L 

•  I ) *DAR ( I  ) 

..  _ -X2  .CONTINUE . .  __  _ _ _ 

SP1  =  P<1)*P(3>*2.*P(6)*2.*PCB>-RHN*<2.*P<10>*P<4)> 

_ S £2 _ »  P(2)*P(3)*P(4I*P(5)*2.»P(7>*P(8)*2.»P(9)-R0N»<2.»P(1Q)*P(4 


*) ) 


_  _S£3  _?  0 , . . 

DO  13  I ali 10 
U_SR3l  a_SP3  __*p  ( I) .. .  . 
SP3aSP3"PR*PAR 


_ c _ _ _ 

C  SOLUTION  OF  MATRIX 

_ P.E.N?_  .  4AJD*  l.A  2  (2  )*A3jl_3J."  A3  (2)*A2(3))-A2(1)*(A1(2)»A3(3)-A3(2)» 

*A1(3))*A3(1)»< A1<2)*A2(3)-A2(2)*A1(3)) 

_  DEN_b1.,/DEN  _  .  __  ....  _ _ _ _  _ _ 

XN1  aSPl*(A3(2)*A2(3)-A2(2)*A3(3))*SP2*(Al(2)*A3(3)-Al<3)*A3(2)  )* 

_ •S.P.3*.(M13.L*_A  2 1 2.L-.  All  21  fMJJJJt _ 

XN2  eSPl*(A2(l)*A3(3)-A3(l)*A2(3) ) *SP2* < A3 <1)» A1 < 3 )-Al < 1 > *A3 (3)  >* 

. . . * S.E3*. ( A  111  ).*  A2.t 3  J_-  A 2.(  1.)  *  A 1  ( 3  >_1 _ 

XN3  =SP1*(A3(1)*A2<2)-A2(1)*A3<2))*SP2*(A1(1>*A3(2)-A3(1)*A1(2)  )* 

_ __»5  P3  •AA2A1)*. A 112  >_-  Al(l )_»  A  2  ( 2 JJ _ ; _ 

C 

_ C _ _ CALCULATION  OF  PARTIAL  PRESS  CORRECTION  AND  NEW  GUESSED  PART  PRESS 

DELl*  XN1  *DEN 


_ -D£L2a_X_M2_tDEN _ 

DEL3aXN3»DEN 

_ _  .PLLLsm.)±HELJ _ 

P ( B ) »P ( 8 ) *DEL2 

_ P  (10).=PX1Q )  *DEL3. 
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SUBROUTINE  PROP  (Concluded) 


C 

7 

FIND  ENTHALPY  OR  ENTROPV  FOR  ITERATION  TO  FIND  TEMP. 

CONTINUE 

17 

I F < IT-1 >18, 6  .17 

CAl L  THERMIP.T.H. CP. F.PR.PAR.HTM. 2» 

ia 

GO  TO  00 

C1IL  THERM  IP.T.H.rP.F.PR.PAR.WTM.ni 

80 

CONTINUE 

IF  (1CT.GT.O)  00  TO  85 _ 

DELT*(HO-H)/CP 

_ IF  (DEL  TpDEL.TI  1  fll  .81. M2. _ 

82 

T»T*DELT*,83 

GO  Ttl  R5 

81 

T*T*DELT*.30 

ICT*2 

■ 

85 

1CT*1CT«1 

dpi  Ti  «nEi  t _ , _ 

IF  (T.LE.300.)  T>300. 

_ if  it. be. soon,  i  T»5non. _ 

C 

C 

ITERATION  CQNTRDL--TFMP. WITHIN  1  DEB  OP  PRESS  WITHIN  QEU 

6 

IF ( ABS(HO-H)  -  CP)6.6»4 

21 

22 

I F ( ABS( DELI )  *  DELL)  22*22.4 

IF  I ABSI DFL2I  -  DF1L1  23.23.4 

23 

4 

IF (ABSIDEL3)  -  DELL)  5.5,4 

CONTINUE 

5 

CONTINUE 

PHOLD«PR 

C 

S*H 

CALL  THERM(P,T,H,CP.F,PR,PAR.WTH,2) 

S=H*. 999331 

C 

CALCULATE  ENTHALPY  AND  ENTROPY  ETC.  FOR  TEMP  KN0WN--IT»1 

CALL  THERM (P. T.H.CP.F.PR.PAR.WTM.O) 

H>h*. 999331*9. /5. 

_ Cfi»CP_*.  999131 _ 

WTM*WTM 

_ _ 

NT  aKOUNT 

_ GALL._TJ_HJEfl...<  TLI2J _ 

TIME  ■  TT2-TT1 

_ 8EBISJJ _ 

99  FORMAT  <  41H  ITERATION  DID  NOT  CONVERGE  IN  50  STEPS  ) 

97  FORMAT ( 36H  THIS  CASE  HAS  TOO  LITTLE  0Z  TO  RUN  ) _ 

END 
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SUBROUTINE  THERM 


SUBROUTINE  THERM  < ALP»T,H,CP.F»PR,PAR.WTM,KT> 

_ QLMEN.S LON _ ALP<  10 )  .CCm. CT(14Q),TT(5)  .F(10)- _ 

DATA  CT  /  2. 5»4*0.»2. 5470497 E4*«4» 60 Q1096B-1* 

1*5+4 •  n..g.547n497E4.-4.6  0niC96E-P3.02l8B94,-2.173?249E-3J 


3.7542203E-6,-2,9947200E-9,9.0777547E-13,2.9137190E4(2.«46007«EO, 
2. 5372S67F0.-1 .84921 906-5.-8. 8Q17921E-9. 3. 9643621E-12^ 


5.574360BE-16.2.9230007E4.4.9467942EO,3.B23470BEO»-1.1187229E-3, 
1 . 24666 1 9 E.-.  6j  jl2.«.1  0.3  5A9. 66-1 0_,_^5j.23A6.5  51E^14j.  3.  ..5  S3  2  ZB  7£3j_ 


5. 8253029E-1»2. 8S95544E0 ,9 .9835061E-4, -2 .1B79904E-7, 1 . 9802785E-11. 
-.3 . 8  4  52 9 40.EM. 6. .3., .88 11  742  E.3  l5  *5  547.01 6  Efl_.A,lA6  9 4  7  6Ei)  t  -  j,  1197237E;iJL^ 
9.6922467E-6#-7,B633639E-9»2.23095l2E-12»9,7447894E3,2.5694290EO, 

3 .1529360 EO .1 .40.5 9955E- 3.  -5  .7fl7B462E-7 . 1 . 062B209E-10 . _ 

-7.3720783E-15,9,8522048E3,6.9446465EO,3.7871332EO.-2.l709526E-3i 
S+0 75 73 3 7 E.-j6 , -  3j_4.7_ 17 726E-9,  7.  J2 1 68,41  E - 1 3 , - 1 . 4 3 6350BE4, 2_.6334459E0, . 


2.  95115 19 E0.1. 5525567E-3-, -  6, 1911411E-7, 1, 1350336E-10, 
-7. 78B2732E-15. -1.423 18 2 l£A*4+3  3 14450E0^2 .« 


-9. 6119332E-6, 9. S122662E-9, -3, 3093421E-12, -9. 672537262.^^^1^^^7^ 


*2.4593791E-15,-B.5491002E2,-1.6481339EO»3.7189946EO» 
l-_2 .516728  8E.-3.X  B_,3 8  3 73531- 6^  -.6 .2498716E-9..2  ^7_Qfl2180E-12. 


-1. 0576 706E3, 3 . 9080704E0 . 3.5976129E0 1 7 . 8145603E-4,-2 . 238667 0E- 7 • 
4^2490l59E -11^-3. 34602 0 AE-15 ,  - 1 .1 9 279 1 8 E3  ^3 .7492639E1L. _ , _ 


4. 15650 16 E0,»1. 724433 4E-3, 5. 6982316E-6,-4.5930044E-9, 1.42 3365 4E-12 

i26fl7-70E4j.-6j861^24^E-li2.6707?32EQi31 031711^E^L _ 

-8.5351570E-7,1.1790853E-10#-6.197356BE-15#-2.9B88994E4, 

6 . 8  838 341 E  0  .  17  0  IQ  QJEJjJ.  .0  3 781  15E-2trJUP7jL39  38  E -5,  6..3_4  5 9 1 75  E - 9 . _ 

-1.628o7oiE-12f -4.8352602E4»1.066438BEl»4.4129266E0#3.l922B96E-3f 

■1.29782  3  0  E-  6j2-»A1  47448  E-10  ,_-P6  74 2 9  B  6  E  - 14  ,-4_,8944Q43E4 , _ 

•7, 2875769E-1, 3. 69 16148E0, -1.333 2552 E- 3,2.650 3100E-6, 

■9 . 7688341E- 10 . -9 . 9772234E- 14 . - 1 , 0628336E3. 2 . 2874980E0 . 2 . 8545761 


,1.5976316E-3,-6.2566254E-7.1.l315849E-10,-7.6897070E-l5, 
-8j_9  0 1 7445 E2 , 6 . 39  Q2879E0/ 

TTLN  «  ALOQ(T) 

TLN  ■  l.-TTLN _ 


CP>0. 


_ Hiflj _ 

RT«1 ,9871B*T 

TTfl)aT _ 

DO  1  K»2,4 

_1_  TUJ4  )_«_T_TJK  rll*T _ 

TTI5)*1,/T 

_ ka*z _ : _ 

IFCT.LE.1000,)  KA*0 

_ DD__2lJS»J,-1D _ 

KK»  14*  (K-P+KA 

_ _ 

KKK«(KK*J) 

3  CC( CT (KKK  )4TTtJ-H _ 

CC ( 1 ) a  CT(KKfl) 

_ LE(KT-l)  10  <  4 1.5 _ 

10  F ( K  )■  CC(1)*.5*CC(2>*CC(3)/3.*.25*CC(4)*.2*CC(5)*CC(6) 

_ 0.Q_.IQ_8 _ _ 

4  F ( K ) ■ ( CC ( 1) *TLN  ) - . 5*CC (2) - <CC< 3 >♦ .5*CC( 4) ) /6, - . 05*CC(5 )*CC( 6)- 

_ «£T_CKK*7J _ 

QO  TO  2 
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SUBROUTINE  THERM  (Concluded) 


5  IF  (ALP(K) .LE.O. )  BO  TO  2 

_ E.CK CC  { 1 ) •TTLN^CCi 2J  ♦  t3«.CC.(l) ♦C.CiiJ _ 

• .  25*CC  (5  >  *•  CT  <KK+7)  -  ALO(3(  ALP  (  K  )  ) 

H  ALP(K)*F(K1 _ _ 

CPP  >  ALP (K )* (CC( 1  )+CC <  2  )+CC ( 3 ) *CC ( 4  )  *CC ( 3 ) ) 

_ CP.  *  .CP  *..CPR - - - - - - 

2  CONTINUE 

_ PAHL»-0*. _ _ _ _ _ 

IF  (PAR.GT.O.)  PARL«ALOQ<PAR) 

_ MTit  «  ALP_<1)«1. 00797  *ALP< 2 )*15 . 9994»ALP< 3>*17 . 0074*ALP< 4 ) *30 . 0061 

1*ALP  <5)*28»01r'55*ALP<6)*2.0l594*ALP<7)*31.998B*ALP(O>*lfl*0l534 

,.2:*ALP.{91*44. 00995*ALP.(10)*28  •  0134*fiAR*39..94B_ _ _ _ 

CP=(CP*PAR*2.5)*<1.98718/HTM) 

_  LF(KT-l)  7,8.9  _  . . . . . 

7  Hs<PAR*<2.5-<745.375/T))*H>*<RT/WTM) 

_ _ _ _ _ 

QO  TO  8 

J2_Sj*_LH_*.  PAR <ti 2.5 •  TJUN_*4jJ6 6107 6  ,-PARl  ))ttl.9B718yMIMl _ 

DSOT  aCP/T 

UsC 

“CP-DSOT  '  "  "  . .  “ 

.  _WTM«MTM/PR  _  _ _  .  _  _ _  _ _ 

8  RETURN 

_ END _ _ 
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COMPUTER  PROGRAM  CALCULATIONS 
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TABLE  I  (Continued) 


specif*  enthalpy  equivalence  ratio 

AND  PRESSURE 

fOR  A  Ch2 

HY0ROCAR8ON 

NITh  air 

.Ho.  .  . 

P 

T 

H 

s 

CP 

HOL 

NO 

TYPE 

TIRE 

ratio 

ATN 

DEG  R 

9TU/L0 

BTU/L0  OR 

BTU/LB 

OR  WT  ■ 

ITER 

CASE 

SEC 

1.1000 

.'i 

"4073.1 

.37”" 

47555  " 

.34626 

20.058 

20 

0 

3.425 

1.1000 

•  i 

4323.5 

200.11 

2.52307 

.34705 

27.000 

21 

0 

3.506 

1.1000 

■  i 

4*2? .  i 

400-08 

—r.'5ffl'2'4 

•  38107 

ST.  255 

Si 

0 

T.50T 

1.1000 

.i 

4091.3 

600.21 

2.61166 

.34900 

20.815 

23 

0 

3.915 

1.1000 

.1 

4038.7 

799.78 

2.65354 

.35039 

26.352 

29 

0 

4.900 

1*1000 

1*0 

4206.4 

.41 

2.31307 

.34759 

20.214 

12 

0 

2.113 

1.1000 

1.0" 

4536.4 

200.05 

2.35067 

.34963 

27.940 

21 

0 

3.509 

l.iono 

1.0 

4793.6 

400.05 

2.40152 

.35077 

27.587 

21 

0 

3.507 

1.1000 

1.0'  " 

$011.5 

600.19 

2.44234 

. 35156 

27.194 

13 

0 

2,274 

1.1000 

1.0 

5202.9 

799.09 

2.40143 

.35220 

26 . 779 

12 

0 

2.110 

1.1000 

10.0 

4203.0 

.37 

2,15120 

.34B34 

28,302 

10 

0 

1.704 

1.1000 

10.0 

4707.7 

199,95 

2^19561 

.35104 

28.141 

15 

0 

2.604 

1.1000 

10.0 

5037,3 

400.00 

2.21664 

.35248" 

27.874' 

21 

0 

3.509 

.  1.1000 

10.0 

5315.5 

600.13 

2.27529 

.35330 

27.552 

14 

0 

2.437 

1.1000 

io.o 

$560  .“4 

800.00 

2.31203 

.35403 

27.199 

17 

0 

2.929 

1.5000 

,1 

3737.7 

.30 

2.57079 

.35452 

26.343 

7 

0 

1.293 

1.9000 

.1  ■ 

4170.9 

199.88 

2.62111 

.35699 

26-213 

11 

0 

1.950 

1.5000 

.1 

4402.5 

399.96 

2.6672? 

.36154 

25.938 

15 

0. 

2,606 

1.5000 

.1 

4691.3 

600.31 

2.71085 

.36309 

25.568 

12 

0 

2.112 

1.5000 

.1 

4052.3 

000.32 

2.75274 

.36423 

25.157 

20 

0 

3.423 

1.5000 

1.0 

3756.6 

™  T31 

2.39730 

.85 vsr 

““  2578*6 

9 

0 

1.622 

1 .5000 

1.0 

4264.7 

199.79 

2.44702 

.35959 

26.313 

7 

0 

1.291 

1.5000 

i.o 

'  '4677.2 

400.01 

•2.49173 

.36274 

26.159 

12 

0 

2.113 

1.5000 

1.0 

4979.3 

600.36 

2.53317 

.36463 

25.891 

16 

0 

'  3.090 

1 . 50n0 

“'5207.6  " 

"  800.34 

2.57240 

.36500 

25.551 

10 

0 

1.702 

1.5000 

10.0 

3763.1 

.30 

2.22390 

.35473 

26.373 

11 

0 

1.954 

1.5000 

10.0 

4301.1 

“  199.00 

”77273*5 

.35984 

26.399 

11 

0 

1.99Z 

1.5000 

10.0 

4793.4 

400.02  • 

2.31742 

.30342 

26.290 

8 

0 

1.456 

1.5000 

10.0 

5202.9 

600.40 

2. 35747 

.36500 

26.141 

12 

”<r 

2.113 

1.5000 

10.0 

5525.9 

000.39 

2.39472 

.36735 

25.904 

18 

0 

3.097 

2.0000 

.i  ‘ 

3i7S. 9 

759 

”2.6 61$ 6  ~ 

. 36iB6 

24.369 

9 

0 

r  •  827 

2.0000 

.1 

3607.5 

199.81 

2.71973 

,36067 

24.340 

7 

0 

1.293 

2.0000 

.1 

4124.7 

399.91 

2.77090 

.87373 

24.233 

10 

0 

1.705 

2.0000 

•  l 

4456.2 

600.49 

2.01750 

.37699 

24.019 

12 

0 

2.112 

2.0000 

.1 

4699.3 

800.43 

2.06120 

.37925 

23,726 

11 

0 

1.947 

2.0000 

1.0.  . 

3176.2 

.39 

2.47394 

.36109 

24.371 

11 

0 

1.953 

2.0000 

1.0 

37o7'.l 

‘  t99".  0 1 

2.53190 

.36904 

24.361 

9 

0 

1.623 

2.0000 

1.0. 

4201.7 

399.91 

2.58257 

.37423 

24.317 

8 

0 

1.457 

2.0000 

1.0 

4626.0 

600.45 

— 276  7797 — 

, 37790 

24.202 

¥ 

9 

T7TT9 

2.0000 

1.0 

4962.7 

000.56 

2.06906 

.30047 

24,003 

12 

0 

2.112 

2.0000 

10.0 

3176.9 

.39 

2.20632 

.36190 

24.372 

12 

0 

2.120 

2.0000 

10.0 

3714.0 

199.01 

2.34430 

.36910 

24.369 

10 

0 

1.790 

27ffoofl 

io.o 

'4234,6 

’  399.92 

2.39470 

'  .37444 

24.353 

10 

0 

1  •  70B 

2.0000 

10.0  _ 

4720.7 

600.46 

2.43949 

.37639 

24.305 

9 

0 

1.022 

2.0000 

10.0 

5149.3 

600.99 

2.46007 

,38130 

24.'2C2 

— 8 

■  0 

1  i  95© 
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TABLE  I  (Continued) 


SPECIFY  ENTHALPY  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  A  CM2  HYDROCARBON  WITH  AIR 


EQ 

p 

T 

SPECIES  MOL 

FRACTIONS 

RATIO 

ATM 

DEO  R 

H 

0 

OH 

NO 

CO 

.500 

.10 

2802.4 

.1743- 

08 

.8520-05 

.1578-03 

.1162-02 

.7481-05 

.500 

.10 

3407.5 

.2148- 

U« 

.2777-03 

.1451-02 

.3984-02 

.3460-03 

.500 

.10 

3913.0 

.3824- 

03 

.2228-02 

.5341-02 

.8236-02 

.3265-02 

.500 

.10 

4274.7 

■1»33- 

02 

.7304-02 

.1087-01 

.1236-01 

.1058-01 

.500 

.10 

4538.7 

.5174- 

02 

.1532-01 

.1633-01 

.1581-01 

.1978-01 

.500 

1.00 

28n3.1 

.3116- 

07 

.2705-05 

.8900-04 

.1163-02 

.2376-05 

.500 

1.00 

3418.0 

.4101- 

05 

.9233-04 

.8446-03 

.4057-02 

.1159-03 

.500 

1.00 

3973. B 

.9293- 

04 

.8740-03 

.3486-02  _ 

.8894-02 

.1349-02 

.500 

1.00 

4423.9 

.6481- 

03 

.3555-02 

.6259-02 

.1438-01 

.5874-02 

.500 

1.00 

4770.3 

■2112: 

07 

.8695-02 

.1390-01 

.1941-01 

.1363-01 

.500 

10.00 

2803.5 

.5559- 

08 

.8574-06 

.5013-04 

.1184-02 

.7531-08 

.500 

10.00 

3422.6 

.7529- 

08 

.2986-04 

.4825-03 

.4091-02 

.3757-04 

.500 

10.00 

4003.8 

.1*22- 

04 

.3067-03 

.2111-02 

.9234-02 

.4840-03 

.500 

10.00 

4518.9 

•l69i- 

03 

.1456-02 

.5579-02 

.1575-01 

.2612-02 

.500 

10.00 

4950.4 

.7309- 

03 

.4173-02 

.1054-01 

.2242-01 

.7691-02 

.900 

.10 

3927.9 

.7816- 

03 

.1146-02 

.5113-02 

.4035-02 

.1169-01 

.900 

■  10 

4214.9 

.2703- 

OS 

.3315-02 

75582-02 

.6238-02 

.2547-01 

.900 

.10 

4438.6 

.6217- 

02 

.6997-02 

.1450-01 

.8495-02 

.3960-01 

.900 

.10 

4624.2 

.1144- 

01 

.1221-01 

,1936-01 

.1062-01 

.5196-01 

.900 

.10 

4785.0 

.1841- 

01 

.1890-01 

.2377-01 

.1254-01 

.6205-01 

.900 

1.00 

4034.3 

.2444- 

03 

.4946-03 

,3578-02 

.4328-82 

.6260-02 

.900 

1.00 

4396.5 

.1092- 

JL2 _ 

.1665-02 

.7455-02 

.6963-02 

.1735-01 

.900 

1.00 

4882.1 

.2924- 

02 

.3914-02 

.1222-01 

.9916-02 

.3086-01 

.900 

1.00 

4920.4 

.5978- 

02 

.7412-02 

.1735-01 

.1286-01 

.4400-01 

.900 

1.00 

5128.1 

.1046- 

01 

.1226-01 

.2245-01 

.1567-01 

,5559-01 

.900 

10.00 

4100.4 

.6056- 

JL2 _ 

.1861-03 

.2277-02 

.4478-02 

.2711-02 

.900 

10.00 

4?40.7 

.35/7- 

o3 

.7208-03 

.5258-02 

. 74oi-bi 

.9922-02 

•  90Q_ 

10.00 

4897.5 

.1152- 

02 

.1889-02 

.9340-02 

.1085-01 

.2114-01 

.900 

10.00 

5198.6 

.2859- 

02 

.3897-02 

.1414-01 

,1456-01 

.3375^51 

.900 

10.00 

5480.1 

.5088- 

02 

.6911-02 

,1933-01 

.1832-01 

.4604-01 

1.000 

.10 

4035.8 

.1*13- 

02 

.1187-02 

.5657-02 

.3251-02 

.2436-01 

1.000 

.10 

4286.0 

.4214- 

M _ 

.3299-02 

.1012-01 

.5345-02 

.3906-01 

1.000 

.10 

4490.6 

.8497- 

02 

.6828-02 

. 1501-01 

.7486-02 

.5281-01 

1.000 

.10 

4864.1 

.1450- 

01 

.1174-01 

.1976-01 

,9488-02 

.6435-01 

1.000 

.10 

4816.5 

.2225- 

01 

.1801-01 

.2404-01 

.tiio-oi 

.7362-01 

1.000 

1.00 

4181.8 

.8269- 

03 

.4909-03 

.3944-02 

.3133-02 

.1714-01 

1.000 

1.00 

4494.2 

.1903- 

02 

.1619-02 

77879-02 

.5662-02 

.3053-01 

1.000 

1.00 

4753.3 

.4282- 

02 

.3780-02 

.1266-01 

.8476-02 

,4448-01 

1-000 

1.00 

4975.0 

.7959- 

02 

. 7102-0i 

.1778-01 

.1127-01 

.5718-01 

1.000 

1.00 

5172.1 

.1313- 

01 

.1187-01 

.2281-01 

.1393-01 

.6803-01 

1.000 

10/00 

4294.8 

.2088- 

03 

,1685-oT 

.2450-02 

.2721-02 

,ioW-oi 

l.ono 

10.00 

4872.9 

.7366- 

03 

.6629-03 

.3476-02 

.5427-02 

.2176-01 

1.000 

10.00 

4993.2 

.1861- 

02 

.1769-02 

.  9631*02 

TB  762-02 

.3456-01 - 

1.000 

10.00 

5271.3 

.3769- 

02 

.3659-02 

.1444-01 

.1231-01 

.4735-01 

1.000 

10.00 

5521.1 

.6692- 

OS 

.6500-02 

""  .  1943-01 

.1590-01 

.5926-01 
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TABLE  I  (Continued) 


Specify  enthalpy  equivalence  ratio 

n7Ii^J:I4-4cITI:iaa.]:WWRTP»n7n:7iT^PTi1Il  111  ■  II  —Mi— 

Ed 

P 

T 

SPECIES  H0L 

fractions 

v 

RATIO 

ATH 

DEO  ff 

H 

0 

OH 

"  NO 

CO 

l.mo 

.10 

4073.1 

.2475-02 

.8478^03“ 

.4977-02“ 

.2123-02 

.4143^01 

1.100 

.10 

4323.5 

.5787-02 

.2840-02 

.9752-02 

.4231-02 

.5442-01 

1.100 

.10 

4522.7 

. 1082-01 

.6160-02 

■man 

MINIUM 

1-100 

.10 

4691 • 3 

.1763-01 

.1060-01 

.1952-01 

,8266-02 

.7705-01 

i.ioo 

.10 

4836.7 

.2610-01 

. 1674-01 

7237  5-01 

.1000-01 

.8525-01 

1.100 

1.00 

4706.4 

.1040-02 

.2535-03 

.2934-02 

.1522-02 

.3619-01 

l.mo 

1.00 

4536.4 

.2761-02 

.1244-02 

.7187-02 

.4005-02 

.4714-01 

1.100 

1.00 

4793.6 

.5681-02 

.3246-02 

.1220-01 

.6782-02 

.5954-01 

1.100 

1.00 

5011.5 

■  II  1  M  1  11 

.6375-02 

.1742-01 

.9527-02 

.  .7102-01  “ 

1.100 

1.00 

5202 .9 

.1584-01 

.1064-01 

.■  .2245-01 

.1209-01 

.6071-01 

1.100 

10.00 

4263.0 

.3796-03 

.4644-04 

.1308-02 

,7605-03 

.3307-01 

1.100 

10.00 

4707.7 

.1150-02 

.3988-03 

.4409-02 

.3059-02 

.4013-01 

1.100 

10.00 

5037.3 

.2603-02 

.l366-0T~ 

.8798-02 

.6309-02 

.5096-01 

1.100 

10.00 

5315.5 

.4942-02 

'  .3112-02 

.1382-01 

.9819-02 

.6248-01 

1.100 

10.00 

tUXKl 

—mime 

1.500 

.10 

3737.7 

.1959-02 

.1303-04 

.5090-03 

.6863-04 

.1189*00 

1-500 

.10 

4178.9 

.7575-02 

.3066-03 

.3032-02 

.5882-03 

.1216*00 

1.900 

.10 

4482.5 

.1665-01 

.1759-02 

.7924-02 

.1878-02 

.1245*00 

1.500 

.10 

4691.3 

.2742-01 

.4706-02 

.l323^0T 

.3479-02 

.1279*00 

1.500 

.10 

4852.3 

.3932-01 

.8844-02 

.1791-01 

.5026-02 

.1311*00 

1.500" 

3756.6 

.6632*0^ 

.1511-05 

.1756-03 

■ *402-04 

.1192*00 

1.500 

1.00 

4264.7 

.3057-02 

.5345-04 

.1331-02 

.2730-03 

.1222*00 

1.500 

4677.2 

.8342-02“ 

•  52B3*  0 3 

,4805*02 

.1279-02 

.1247*00 

1.500 

1.00 

4979.3 

.1986-01 

.2049-02 

.1003-01 

.3094-02 

.1278*00 

1-500 

1.00 

5207. 8 

■ni  r  tjim 

cssgiTT  r/1  CTTH1 

■ft™ 

Ml  II  1  1  1  I  I 

•  l3l2*Tfo 

1*500 

10-00 

3763.1 

.2147-03 

. 159Q-Q6 

.5722-04 

.7861-05 

.1193*00 

1.500 

10.00 

mm 

.1069-02 

Ml  II  1  1 1 1  ■ 

Ml  1  1  1  1  1  1  1 

.  100*-  D<S 

l<5no 

10*00 

4793.4 

. 3405«02 

.9704-04 

.2181-02 

.6189-03 

•1248*00 

1.500 

10-00 

5202.9 

“  .7633-02 

• 2020“ 0 ^ 

.5*10-0* 

-  -.2050-02“ 

•  1273*00 

1.500 

10.00 

5525.9 

.1345-01 

.1848-02 

.1104-01 

.4312-02 

•  1303-00 

2.000 

.10 

3173.9 

.2984-03 

. 1849-07 

.1071-04 

.5017-06 

7T791*00 

2.000 

.10 

3687.5 

.2519-02 

.2476-05 

.1767-03 

.1407-04 

.1826*00 

.10 

.1012-01 

^1  1  1  1  1  1  1  ■ 

.10*0-0* 

—  1  1  1  1  1  1  1  I 

—  1  1  1  1  1  1  1 

2.000 

.10 

4456 . 2 

.2396-01 

.4615-03 

.3372-02 

.4878-03 

.1837*00 

2.  000 

.10 

“4OT9.3 

.4113-01 

. 1634-02 

. 0056-9* 

.1133-02' 

.1826*55 

2.000 

1.00 

3176.2 

.9541-04 

.1894-08 

.3433-05 

.1613-06 

.1791*00 

2.000 

1.00 

"3707.1 

t  9957*0^ 

72896- oS- 

•  6122*04 

- 74947-05 

.1831*00 

2.000 

1.00 

4201.7 

.4024-02 

.1000-04 

.4611-03 

.5487-04 

.1852*00 

2.000 

l.Ao 

MU  1  1  1  1  1  ■ 

,1136-03 

.1810-0* 

, *845-03 

^■i  mini 

2.  000 

1.00 

4962.7 

.2351-01 

.5795-03 

.4473-02 

.8499-03 

.1854*00 

2.000 

10.00 

3176.  r~ 

. 302Bs04 

.1909-09 

.1091-05 

.5127=TTT 

•  1792 400 

2.000 

10.00 

3714.0 

.2775-03 

.3058-07 

.1998-04 

.1623-05 

.1832*00 

TABLE 


(Continued) 


SPECIFY  ENTHALPY  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  A  CH2  HYDROCARBON  KITH  AIR 


Ed 

p 

T 

SPECIES  MOL 

FRACTIONS 

RATIO 

atm 

DEG  R 

H2 

02 

H20 

C02 

N2 

A 

•  900 

.10 

2602. * 

.2*38-05 

.1006*00 

.6739-01 

.6779-01 

,7939*00 

.9009-02 

•  900 

.10 

3407.9 

I^H  T'TT  D'LBi 

.9887-01 

.6661-01 

.6737-01 

.7920*00 

.9003-02 

.900 

.10 

3913.0 

.6272*03 

.9621-01 

.6418-01 

.900 

.10 

4274.7 

.1892*02 

.9396-01 

.9823-01 

.9623-01 

.7377*00 

.8883-02 

.900 

.10 

4938.7 

.3543-02 

.9179-01 

.5134-01 

.4616-01 

.8764-02 

.900 

1.00 

2803.1 

.8378*06 

.1006*00 

.6743-01 

.6776-01 

.7539*00 

.9009-02 

.900 

1.00 

3418.0 

.2741-04 

.9897-01 

.6700*01 

.9006-02 

_ -soo 

1-00 

3973.8 

.2520-03 

.9606-01 

.6529-01 

.6627-01 

.7484*00. 

.8990-02 

UHifiia 

1.00 

■  1  11  IT  — 

n/imi 

1MI-PTW 

—  1 II 1  II  !■ 

■uim  i  — 

■LUIOMIH 

.900 

1.00 

4770.3 

.2209-02 

.9039-01 

-  .9604-01 

.9299-01 

.7317*00 

.8853-02 

.900 

10.00 

2803.9 

.2656-06 

.1006*00 

.6745-01 

.6776-0! 

.7539*00 

.9009-02 

.900 

10.00 

3422.6 

.8882-09 

.9904-01 

.6721-01 

.6772-01 

.7524*00 

.9008-02 

.900 

10.00 

4003.6 

.8955-04 

.9611-01 

.6629-01 

.6721-01 

.7492*00 

.9001-02 

.900 

10.00 

4918.9 

■rttb™ 

■WH-II 

■ .6393*01 

.6490-01 

.7438*00 

.8976-02 

.900 

10.00 

4950.4 

Kf  if  ')>  i 

■frl-VJiOTM 

■mini  — 

■MfTTPT 

.900 

.10 

3927.9 

,2382-02 

.2289-01 

.1118*00 

.1057*00 

.7257*00 

.8690-02 

.900 

.10 

4214.9 

.5162-02 

.2811-01 

.1044*00 

,90ffi-Dl 

.7159*06 

.6585-02" 

.900 

.10 

4438.6 

.8390-02 

.3316-01 

.9932-01 

.7471-01 

.7042*00 

.8499-02 

.900 

.10 

4624.2 

.1157-01 

.3707-01 

.8823*01 

.6050=01 

.  6917*00 

.8323=02 

.900 

.10 

4789.0 

.1457-01 

.3963-01 

•  .7462-01 

.4849-01 

.6788*00 

.8180-02 

.900 

1.00 

4034.3 

.1193-02 

.2009-01 

.1117*00 

■Hj-LID.LM 

.900 

1.00 

4396.9 

.  .3178-02 

.2379-01 

.1092*00 

,9996-01 

.7211*00 

.8652-02 

.900 

1.00 

4682.1 

inmM 

■UiUW 

,8466-01 

.7110*00 

.900 

1.00 

4920.4 

.8717-02 

.3232-01 

.9328*01 

.6994-01 

.6997*00 

.6432-02 

.900 

1.00 

■94»i-IJ<« 

■  1  IIHIIM 

,9669-01 

.6877*00 

.8309-02 

.900 

10.00 

4100.4 

.4964-03 

.1838*01 

.1163*00 

.1196*00 

.7308*00 

.8752-02 

10.00 

■  H  N  1  1  1  ■ 

.2007-01  ‘ 

■U11Q1H 

■HMfllW 

.900 

10.00 

4897.9 

.3553*02 

.2332-01 

.1079*00 

.9549-01 

.7172*00 

■  8628-JI2 

.900 

10.00 

5196.6 

.9099-02 

.2690-01 

,1007*90 

■W'ttP'IM 

.900 

10.00 

5460.1 

.8993-02 

.3006-01 

.9284*01 

.6784-01 

.6966*00 

.6427-02 

,  1.000 

.10 

4035.6 

W-HTHM 

.1194*00 

.1041*00 

1.000 

.10 

4286.0 

.8497-02 

.1798*01 

.1108*00 

.8770-01 

.7045*00 

.8444-02 

1.000 

.10 

4490.6 

■  H  1  1  1  1  !■ 

wr+t  DM 

■  1  1  1  1  1  1  1 

1.000 

.10 

4664,1 

.1851-01 

.2800-01 

.8986*01 

.9841-01 

.6802*00 

.8176-02 

1.000 

.16 

"  "4816. 5 

.1169*81' 

“T3103-01 

.7689-01 

,4703-01 

.6674*00 

,8037-02 

1.000 

1.00 

4181.8 

.3345-02 

.7600-02 

.1234*00 

.1121*00 

.7196*00 

.8611-02 

1.000 

DUJ 

,1168*00 

hi  I  M ii 

1.000 

1.00 

4753.3 

.9063-02 

.1854-01 

.1084*00 

.8176-01 

.7001*00 

.  .6410-02 

— 1 ^■n  1 

■fill 

KDi-ITMi 

■  Mill  1  !■ 

HU 

1.000 

1.00 

9172.1 

.1870-01 

.2664-01 

.8863-01 

.5452-01 

.6768*00 

.8164*02 

1.000 

10.00 

n  rnr 

■  1  1  1  1  1  1  1  — 

.1189*00 

1.000 

10.00 

4672.9 

.3864-02 

.8982-02 

.1216*00 

.1071*00 

.7162*00 

.8584-02 

1.000 

10.00 

4993. Z 

■  1  1  Nil  1  ■ 

■  Mill  — 

■LLLfAUMIMi 

1.000 

10.00 

9271.3 

.9045-02 

.1754-01 

*1076*00 

.7871-01 

.6971*00 

.6398-02 

1.000 

mm 

5521.1 

.1210*01 

.2128-01" 

“  .9891-01  “ 

,65l8-Bl 

.6663*00 

.8290*82 
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TABLE  1  (Concluded) 


SPECIFY  ENTHALPY  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  a  CM2  HYDROCARBON  WITH  AIR 
EQ  P _ T  SPECIES  MOL..  FRACTIONS  4 _ 


RATIO 

ATM 

DEG  R 

H2 

02 

H20 

C02 

N2 

A 

1.100 

.10 

4073.1 

.9758*02 

.4583-02 

.1250*00 

.0736-01 

.7030*00 

,4407-02 

1.100 

.10 

4323 .5 

•1301*01 

.1055-01 

.1159.00 

.6254-01 

.6927*00 

.8296-02 

l'lOO 

■10 

4522.7 

.1665*01 

.1610*01 

■1052.00 

.6017-01 

.6810*00 

.8169-02 

1-100 

•10. 

4691  .3 

.2020*01 

.2057-01 

.9361-01 

.5550-01 

.6687*00 

.6034-02 

1.100 

.10 

4836.7 

.2335-01 

.2383-01 

,8iB2-oi 

.4310-01 

.6563*00 

r7’4’56-02 

1-100 

l.on 

421)6.4 

.7972*02 

.1725-02 

.1293*00 

.1034*00 

.7072*00 

.8453-02 

l.ioo 

1-00 

4536.4 

.1025-01 

.6156-02 

.1227*00 

.*106-01 

•6991*00 

.8372-02 

1.100 

1.00 

4793.6 

•1342*01 _ 

.1127-01 

.1138*00 

.7691-01 

.6669*00 

.8266-02 

1*100 

1-00 

5011-5 

•1692-01 

*1586-01 

.1036*00 

.6349-01 

.6776*00 

.6146-02 

1.100 

1.00 

5202.9 

_ -2041*01 

.1951-01 

.9265-01 

.5179-01 

.6659*00 

.6024-02 

1.100 

10.00 

4283.0 

.7029*02 

.3640-03 

.1310*00 

.1069*00 

.7098*00 

. 8480-02 

1.100 

10.00 

4707.7 

.8053.02 

.2605-02 

.1281*00 

.9906-01 

.7046*00 

.8432-02 

l.mo 

10.00 

5037,3 

, 1036-01 

,6506-02 

.1215*00 

,6692-01 

,6963*00 

.8352-02 

1.100 

10.00 

5SJ5.5 

.1331*01 

.1066-01 

.1133*00 

.7381-01 

.6665*00 

.8255-02 

1.100 

10.00 

5560.4 

.1661*01 

-1438-01 

.1030*00 

.8132-81 

.6790*00 

.8149-02 

1.500 

.10 

3737.7 

.5349*01 

.1238-04 

.1103*00 

.5436-01 

.6448*00 

.7699-02 

1.500 

.10 

417879 

.4977-01 

rsosa-os 

.1171*00 

.5074-61 

.6413*00 

.7661-02 

* 

1.500 

.10 

4482.5 

.4841-01 

.1651-02 

.1096*00 

.4609-01 

.6339*00 

.7581-02 

1.500 

.10 

4691.3 

.4888-01 

.3»05-02 

,4670-01 

.4021-01 

.6241*00 

.7473-02 

1.500 

.10 

4852.3 

.4984-01 

.6288-02 

.0674-01 

.3434-01 

.6132*00 

.7352-02 

1.500 

1.00 

3756.6 

.5373*01 

.1437-05 

.  ii46*ofl~ 

“  75414-01 

.6493*00 

.>706-02 

1.500 

1.00 

4264.7 

.5015-01 

.5369-04 

.1204*00 

.5003-01 

.6439*00 

.7690-02 

1.500 

1.00 

4677.2 

.4817*01 

.5278-03 

.1170*00 

.4732-01 

.6396*00 

.7645-02 

1.500 

1.00 

4979.3 

.4818-01 

.1883-02 

.1089*00 

.4243-01 

.6322*00 

.7567-02 

1.500 

1.00 

5207.8 

.4942-01 

.3810-02 

•  4821-or 

,36lo-0l 

.6226*00 

.7468-02 

1.500 

10.00 

3763.1 

.5381-01 

.1511-06 

.1192*00 

.5413-01 

.6455*00 

.7708-02 

1.500 

10.00 

4301. 1 

.5039-01 

.6770-05 

. 1219*00 

.9001*01 

.(490*05 

.7703-02 

1.500 

10.00 

4793.4 

.4813-01 

.1005-03 

.1217*00 

.4600-01 

.6432*00 

.7684-02 

1.500 

10.00 

5202.« 

.4733-01 

.5844-03 

.1175*80 

.4463-01 

.6388*00 

.7640-02 

1.500 

10-00 

5525.9 

.4809-01 

.1687-02 

.1090*00 

.4001-01 

.6319*00 

.7571-02- 

2.000 

.10 

3173.9 

TITi3»00 

.4753-08 

.7557-01 

.2015-01 

.5767*00 

.6910-02 

2.000 

.10 

3687.5 

.1262*00 

.6696-06 

.7919-01 

.2421-01 

.5760*00 

.6902-02 

2.000 

.10 

“4124.7 

.1204*00 

.1716-04 

.7981-01 

.2202-01 

.9794*00 

.6871-02 

2.000 

.10 

4456.2 

.1133*00 

.1303-03 

.7709-01 

.2056-01 

.5702*00 

.6811*02 

2.000 

.10 

4699.3 

,1061*00 

.4525-03 

.7153-01 

.5624*00 

.6726-02 

2.000 

1.00 

3176.2 

.1314*00 

,4868-09 

.7562.01 

.2012-01 

.5766*00 

.6911*02 

* 

2.000 

1.00 

3707.1 

.1269*00 

.7816-07 

.7998-01 

.2405-01 

f57B9«00 

.6*08-02 

2.000 

1.00 

4201.7 

.1230*00 

.2777-05 

.8136-01 

.2161-01 

.5774*00 

.6895*02 

2.000 

1.00 

4626.0 

.1182*00 

.3196-04 

.0070-01 

.2004-01 

.9746*00 

2.000 

1.00 

4962.7 

.1127*00 

.1616-03 

.7723-01 

.1672-01 

.5696*00 

.6006*02 

2.000 

10.00 

3176.9 

.1314*00 

.4906-10 

, 7963-01 

72811-01 

•9780+00 

.6911*02 

2.000 

10.00 

3714,0 

.1271*00 

.8247-08 

.7971-01 

.2399-01 

.5767*00 

.6910*02 

2.000 

10.00 

4234.6 

.1241*00 

.33*4-06 

.fZoS-TT" 

.2144-01 

.5783*00 

,0*09*02 

2.000 

10.00 

4720.7 

.1210*00 

.5170-05 

.0776-01 

.1977-01 

.5771*00 

.6092*02 

2.000 

10.00 

5149.3 

.1179*00 

.36&z-o» 

.8192-01 

.1892-91 

.9749*00 

.0663**2 

* 

AEDC-TR-71  -256 


Ti 


SPECIFY  TEHPERTURE  EQUIVALENCE  RATIO  and  PRESSUi 
_ EO  P _ T _ H _ 


RATIO 

ATM 

OUTS 

BTU/LB 

.5000 

.1 

3000.0 

63.06 

.5000 

.1 

4000.0 

442.16 

T5ooo 

.1 

_ 

SoOO.O 

1295.10 

.5000 

1.0 

3000.0 

62.49 

,5000 

1.0 

!  4000,0 

4io,5l 

.5000 

1,0 

5000.0 

963,64 

.5000 

HOE 

i  3000.0 

62.21 

.5000 

10.0 

'  4000.0 

398.84 

.5000 

10.0 

m  u  m 

.9000 

.1 

3000.0 

-391.15 

.9000 

.1 

|  4000.0 

*TT5i 

.9000 

.1 

.9000 

1.0 

3000.0 

•391.07 

.9000 

1.0 

_ *16.42 _ 

•  9000 

1.0 

5000.0 

■  N  II  i  M  ■; 

.9000 

3000-0 

•392.32 

.9000 

10.0 

4000,0 

-40.11 

.9000 

10.0 

5000.0 

465.69 

1.0000 

.1 

3000.0 

-498.27 

1.0000 

.1 

4000.0 

•24.84 

1.0000 

9  1 

5000.0 

■  1  II  II  1  — 

1.0000 

1  <0 

3000.6 

-500.95 

1.0000 

1.0 

4600.0 

-98.19 

1,0000 

1.0 

5000.0 

624.31 

1.0000 

10.0 

3000. 0 

-502.21 

1.0000 

10.0. 

4000.0 

-133.00 

1.0000 

10.0 

5000,0 

9  0 j . 11 

1.1000 

.1 

3000.0 

•454.00 

1.1000 

.1 

4000.0 

-47,40 

i.inoo 

.1 

5000. 0 

1047.05 

1.1000 

1.0 

3000.0 

-454.15 

_ 1-10"0  -- 

1.0 

4000 .0 

-96.05 

TTo 

1.1000 

10.0 

3000.0 

-454.20 

1.1000 

10.0 

4000.0 

-108.37 

1.1000 

10.0 

5000.0 

375.55 

1.5000 

.1 

3000.0 

=370786 

1.5000 

.1 

4000.0 

111.49 

1 .5oiO 

.1 

1  bl9 , 7  7 

1.5000 

1.0 

3000.0 

-271.09 

1.5000 

1.0 

4000.0 

92.56 

1.5000 

1.0 

5000.0 

616.42 

1.5000 

1.5000 


II 


RE  FOR  A  CH2 

hydrocarbon 

WITH  AIR 

S 

CP 

H0L 

NO  TYF 

e  time 

BTU/LB  OR 

wT 

SE  sec 

2.28869 

.31243 

28.930 

2  1 

.322 

2.39664 

28 . 754 

2 

.322 

2.58361 

.32648 

'27.663 

3 

•  373 

2.13044 

.31243 

28.932 

2 

.322 

W W*  -M 1 1  •  V 

■n  1 1 1 1  — 

8-rw  TV 

3 

.373 

2,35194 

.32862 

28.141 

3 

.372 

uni  n— 

28.933 

2 

.  322 

2.06887 

,32330 

28.906 

3 

.372 

i  II 1  M  1  — 

'787  525 

3 

.  372 

2.30552 

.32718 

28.906 

3 

.372 

2.42692 

.33890 

28.568 

8  424 

EE3T1  ■ 

.34390 

26.222 

4 

.423 

wnammsm 

umi  m 

2 

.322 

2,25424 

.33925 

28.780 

2 

•  322 

ttn.  1 1 1 

■MUIK 

27 • 63$ 

4 

ram  s 

Kh7\SB 

28.910 

2 

.322 

2.08972 

.33940 

28.862 

3 

.375 

2.20149 

.34521 

28.337 

4 

.423 

2.30733 

.33073 

28.887 

5 

.474 

2.44134 

.34249 

28.427 

5 

.473 

■  min 

26.003 

4 

*4Z& 

2,14817 

.33075 

28.897 

7 

.575 

■  linn  — 

4 

.423 - 

2.42222 

.  .34798 

27.446 

4 

.423 

1.98953 

.33075 

28.901 

7 

.575 - 

2.09516 

.34316 

28.797 

5 

.474 

i  ii  1 1 1 1 

,34890 

28.172 

6 

±  LM&sa&LZ... ... 

2.34839 

,33333 

28.347 

6 

.465 

i  ii 

.34576 

28,132 

“15 

.949 

2.70379 

.35124 

25.771 

7 

.506 

2.1876* 

. 3333S 

ZB  a  440 

4 

'2.28951 

.34591 

28.294 

3 

.335 

1  II  1  1  II  — 

.35152 

"27.218 

6 

a  463 

2.02572 

,33333 

28.348 

4 

.379 

2.i25o3 

.34596 

"287334 

5 

2.23176 

.35233 

27,910 

10 

.634 

^ .  49012 

-.34444 

26 ■ 376 

4 

2.59950 

.35729 

26.288 

4 

,379 

■  iimi  i 

1  1  II  1  1  1 

3 

<42D 

2.31661 

,34444 

26.377 

4 

.379 

2 . 42107 

— 735719 

26 • JiO 

.■379 - 

2.53639 

.36475 

25.866 

5 

.422 

rfl  t  r  r  f  i  ■■kUsit 

,34444 

26.377 

4 

mro 

2.24618 

.35716 

26.869 

4 

.378 

“4733696 

‘  .36*68 

"28  .'231 

3 

•  *22 

> 


D 

9 

9 


K) 

CJI 

0> 


TABLE  II  (Continued) 


infli!  AT 


SPECIES  NOL  FRACTIONS 
- 8 - Oh — 


3000.0 

4000.0 

"TffOo.o 
3000. 0 
'4000.0' 
5000. 0_ 
'Jo'Oo'.o 

4000.0 
"500775“ 
3000 .0_ 
'400  6".  0 
sooo.o 
3000.0 
400 OlO 
Too'oTo" 

3000 -J) 

4000.0 

5000^0 

so'fl'o.o" 

4000. 0 

SOOO.O 

3000.0 

4000.0“ 

5000. 0_ 

3000.0 

4000.0 

5000.0' 

3°0°J.0_1 

4000.0 

5000^0_ 

"3000.0 

4Q00  .0 

"570575“ 
3000  .0_ 
4000.0 
5000. 0_ 
SOOO.'O 
Annn .  n 


1635-0? 

5810-03 

2024-01 

1843-08 

1054-03 

4400-02 

3278-07 

1888-04 

8477-03 

2071-05 

10*0*02 

3218-01 

"35V5-08 

2074-03 

74*6-02“ 

8553-07 

37*5-04 

1558-02 

5502-05 

1392-02 

3535-01 

1195-08 

3030-0? 

8504-02 

2*09-06 

8837-04 

1894-02 

3369- 04 
1*19-02 
3887-01 
1066-04“ 
5511-03 
9734-02“ 

3370- 05 
1703-03 
2384-02 
836i-04 
4543-02 


2646-04 


1853-01 


4588-03 


.3100-04 

.3022-02 
.4484-01 
.9805-05 
.9566-03 
.1487-01 
.3101-05 
.3029-03 
.4*47.-09 
,1361-04 
.1513-02 
.3167-01 
TTSoI-oT" 
.4355-03 
.9053-02 
.1360-05 
■; 1322-03 
.2445-02 
.2125-05 
.1012-02 
.2847-01 
.4509-06 
T»2i-o3 
.7590-02 
.9474-07 
.4750-04 
.1802-02 
.5657-07 
.3611-03 
.2535-01 
.5656-08 
.7038-04 
.6176-02 
.5656-0* 

’  .7437-05 
.1210-02 
"T7B82-0S 
.*352-04 


.7880-0* 


.2222-02 


.9351-0* 


9149-02 

Si76-6i 

1837-02 

9176-02 

2300-01 

1636-03“ 

9194-02 

3394101" 

7959-03 

4!l4-or 

1499-01 

ms=vr 

4125-02 


1205-01 

X23*=*T" 

3003-02 

8315-04 

20*5-02 

1161-01 

5525-04 

1420-02 

8839-02 

3257-05 

1652-02“ 

1184-01 

TO  3  0“ 05“ 

8575-03 

9382-02 

3257-06 


5894-02 
432 3 - 06“ 
2630-03 


1367-06 


3257-02 


2634-04 


4475-02 

3606-01 

*831-05 

1487-02 

2083-01 

3109-0? 

4774-03 

8934.09 

1238-03 

IWFoT 

7328-01 

3*19-04 

5554-02 
4852 -Or 

1240-04 

1695-6*2“ 

2522-01 

8697-03 

2256-01 

8288-01 

4113-03 

11*6-01 

5060-01 

1961-03 

5647-02 

3*89-01 

2616-01 

3854-01’ 

9278-01 

2818-01 

3276-01 

7044-01 

2818-01 

3144-01 

4957-01 

1116*00 

1205*00 


1116*00 


1281*00 


1209*00 
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TABLE  II  (Continued) 


SPECIFY  TEHPERTURE  EQUIVALENCE  RATIO  ANO  PRESSURE  FOR  A  CHS  HYDROCARBON  WITH  AIR 


EO 

P 

_ I _ 

SPECIES ...  HOL 

FRACTIONS 

RATIO 

atm 

M2 

mm 

H20 

C02 

N2 

A 

•  500 

.10 

.9433-05 

aenimn 

.6773-01 

.7535*00 

.6008-02 

.500 

.10 

111  'HI 

■fTTlnTl 

Knifl 

iWiuB'l 

.6287-01 

.7453*00 

.8953-02 

•  500 

.10 

»0V0«Q 

.  8830-01 

.3869-61 

.2518-01 

.6408-62 

•  500 

JrOo.o 

sbf u_rnu 

•1002-OSL 

.6737-01 

.6775-01 

.7536*00 

•9oo9-02 

.500 

1.00 

40  00.0 

.2753-03 

.9569-01 

.6514-01 

■umini 

.8988-02 

5006.JD 

.3432-02 

.8823-01 

.5077-01 

.4508-01 

.7224*00 

.8763-02 

.500 

10.00 

3000.0 

.9450-06 

.1003-00 

.6742-01 

.6776-01 

.7536*00 

.9009-02 

.500 

10.00 

KIiI-f.Hi! 

.8844-04 

.9613-01 

.6626-01 

.6722-01 

.7492*00 

.9001-02 

.500 

10.00 

.1269-02 

.5945-01 

.5851-01 

.8913-02 

.000 

.10 

3000.0 

.3771-04 

.1931-01 

.1180*00 

.1184*00 

.7342*00 

.8771-02 

.900 

.10 

4000.0 

.2948-02 

.2399-01 

•1103*00 

.1026*00 

.7237*00 

.8669-02 

.900 

.10 

5000.0 

.1828-01 

.4111-01 

.5837-01 

.3424-01 

.6589*00 

.7957-02 

.900 

1-00 

3000.0 

.1194-04 

.1929-01 

.1182*00 

.1185*00 

.7343*00 

■mini— 

.900 

1.00 

4noo.o 

•1067-02 

_ .  1988-Ql 

_ .11*6*00 

. 1125*00 

.7293*00 

.8733*02 

.900 

1*00 

Too  o.o 

.9842-02 

r-vl  LJ 1 1  m 

mirw 

Tni.n.f 

SSl£H_  t  Li  U  > 

.900 

1Q.Q0 

3000-0 

.3777-05 

TTTTHU 

.1185*00 

.7343*00 

.8772-02 

.900 

10.00 

4000.6 

.3571-03 

.1168*00 

.1165*00 

.7315*00 

:875B-02 

.900 

10.00 

5000.0 

.4277-02 

.1054*00 

.9097-01 

.7141*00 

.8599-02 

1.000 

.10 

3000.0 

.2662-03 

,4725-03 

.1301*00 

ti-rm™ 

turn.™ 

• 

1.000 

.10 

4000.0 

.4807-02 

.1073-01 

.1203*00 

.1061*00 

.7164*00 . 

.8571-02 

1.000 

.10 

5000.0 

■minw 

nmm 

1.000 

1.00 

3000.0 

.1257-03 

.2121-03 

.1304*00 

.1304*00 

.7297*00 

.8713-02 

1.000 

1.00 

4000.6 

.2290-02 

.5173-02 

ESTWJ  HIM 

.8647-62 

1.000 

1.00 

5000. 0 

.1277-01 

BrttTffT™ 

.6873*00 

.8276-02 

1.000 

10.00 

3000.0' 

.5985-04 

.9362-04 

iTl  l  1 1 1 1 1  ■ 

1.000 

io.oo 

4000.0 

.1092-02 

.2364-02 

.1284*00 

.1247*00 

.7265*00 

.8683-02 

1.000 

10.00 

5000.0 

TtMOU 

Mil '  M  1  1  ■ 

inini— 

1.100 

.10 

3000.0 

.3337-06 

.1299*00 

.1120*00 

.7113*00 

.8493-02 

1.100 

.10 

4000.0 

.9136-02 

.3299-02 

.1268*00 

.1606*00 

.7051*00 

.8429-02 

1.100 

.10 

5000.0 

.2641-01 

.2633-01 

.6750-01 

mrurmm 

■lining 

1.100 

1.00 

3000.6 

rro86-o2 

.1900*00 

mus&mm 

.8494-02 

1.100 

1.00' 

4000. 0 

.7533-02 

.5191-03 

.1312*00 

.1072*00 

.7097*00 

.8477-02 

i.ioo 

1.00 

Mil  1  1  1  1  1  ■ 

Ii»n« 

mm'M 

■  lllllll  — 

1.100 

10.00 

3000.0 

.9987-02 

.3336-08 

,1300*00 

.1120*00 

.7113*00 

.8494-02 

1.100 

io.oo 

4000.0 

.7195-08 

T5798-64 

.1324*00 

. i 087*06 

.71o9«06 

•  8"489-02 

1-100 

10-00 

5000.0 

•1004-01 

.6000-02 

.1224*00 

,8849-01 

.6974*00 

.8362-02 

1.500 

.10 

3000.0 

.6155-01 

.6477-00 

.1116*00 

.6184-01 

.6'45'6*00 

77709-62 

1.500 

.10 

4000.0 

.5118-01 

.9164-04 

.1184*00 

.5238-01 

.6433*00 

.7683-02 

1.500 

.10 

5000.0 

.8547-02 

MIINIIIM 

.2852-01 

■mam 

1.500 

1.00 

.6157-01 

.6475-09 

.1116*00 

.6183-01 

.6456*00 

.7709-02 

1.500 

l.oo 

4000'.  0 

.5191-01 

.9174-05 

.1201*96 

,5252-01 

.6449*90 

, 7701-02 

1.500 

1.00 

5000.0 

.4825-01 

.2023-02 

,1081*00 

.4200-01 

.6315*00 

.7560-02 

1.500 

10.00 

“3600.0"' 

.6158-01 

.6474-10 

"  .1116*09 

.6163-01 

9  6  456^ 0  0 

.7709-62 

1.500 

10.00 

4000.0 

.5217-01 

.9163-06 

.1207*00 

.5254-01 

.6454*00 

.7707-02 

1.500 

10.00 

5000.0 

.4754-01 

.  2578-03' 

.1209*00 

7T659-01 

.6419*00 

; 7666-02 

TABLE  II  (Concluded) 


SPECIFY  TEHPERTURE  EQUIVALENCE  RATIO"  AND  PRESSURE  FOR  A  CH2  HYDROCARBON  WITH  AIR 


EOUL1BRIUN  ISENTRQPJC  EXPANSION  FRON  10  ATH. 


...  EO 

p 

T 

H 

S 

CP 

MOL _ 

NO 

TVPE 

TINE 

RATIO 

ATM 

DES  R 

btu/lb 

BTU/LB  OR 

BTU/LB 

OR  NT 

ITER 

CASE 

SEC 

.5000 

10.(1000 

4000. 0 

398.84 

2.06887 

.32330 

28.906 

3 

1 

.375 

-Sn><0 

S.nooo 

3485.0 

220.84 

2.06887  _ 

.31852 

..  28.927 

7 

.  2 

1  .637 

.5000 

l.-'OOO 

2463.1 

-105.35 

2.06885 

.30335 

28.934 

2 

2 

.569 

.5000 

,5uno 

2104.4 

•213.29 

2.06908 

.29556. 

28, 934 

2 

2 

.568 

.5000 

10. '000 

5000.0 

825.68 

2.16345 

.32903 

28.625 

3 

1 

.374 

.5(100 

S.rooo 

_ 4.4?i.5 

598.16 

2.16345 

.32658 

28.798 

9 

9 

2.064 

.5000 

l.iono 

3319.7 

166.72 

2.16345 

.31663 

28.927 

7 

2 

1.637 

,5000 

,5000 

-  2865.5.  . 

...19.75  . 

_-2. 16344 

,31041. 

28^933 

5 

2 

1.210 

.9000 

10.  ■'000 

4000.0 

-40.11 

2.08972 

.33940 

28.862 

3 

1 

.375 

,9000 

5,0000 

35Q9.6  .. 

-.218.97.. 

.  2.08972.  _ 

33436 

28,900 

7 

7 

1.636 

•  9000 

1.3000 

2515.4 

•549.78 

2.08972 

.31804 

28.910 

5 

2 

1.210 

.9000 

.5000. 

2162.2 

_-66Q,05 

2.08973 

.30961 

28.911 

4 

2 

_ .997 

.9000 

10.0000 

5000. 0 

465.69 

2.20149 

.34521 

28.337 

4 

1 

.425 

. 0000 

5  •  COOO  _ 

.  4565.2. 

_  214.48.. 

.2,20149. 

_ .34319 

28.597, 

_ 12  . 

_ 2  — 

_  2.707  .  . 

.9000 

1.0000 

3512.6 

-214.01 

2.20149 

.33438 

28.889 

7 

2 

1.637 

.9000 

.5000 

3061.1 

•371.35 

2.20149 

.32816 

28.907 

5 

2 

1.210 

1.0000 

10.0000 

4000.0 

•133.00 

2.09516 

.34316 

28.797 

5 

1 

.476 

l.flOOD 

S.roOO 

3536.2. 

-312.83 

_  2  lO«516 

J3B4JL 

28.868 

.2 _ 

1.852 

1.0000 

l.oono 

2556.4 

-648,07 

2.09514 

.32225 

28.904 

4 

2 

.992 

1.0000 

.5000 

2201.6  . 

.-761.17 

2,00515.. 

,31373 

_ 2®  t?Q5 

4 

.  j*04 

1.0000 

10. "000 

5000.0 

405.11 

2.21418 

.34890 

28.172 

6 

1 

.526 

1.0000 

S.r.noo 

4577.7 

172.26 

2.21*18  _ 

.34694 

28j_448 

8 

2 _ 

1.850 

l.oono 

l.-'OOO 

3500.9 

•283.51 

2.21416 

.33898 

28.832 

11 

2 

2.492 

1.0000 

.5000 

3159.5 

•444.45 

2.21416 

.33328 

28.885 

10 

2 

2.276 

1.1000 

io.:oco 

4000.0 

-108.37 

2.12503 

.34596 

28.334 

5 

1 

.423 

1.1000 

5.P000 

3483. 8 

-289.92 

2.12503 

.34036 

28.346 

_ 11 _ _? _ 

2.240 

1.1000 

l.-IOOO 

2488.1 

•623.88 

2.12503 

.32346 

28.348 

11 

2 

2.244 

..  .  .  1.1000 

.5300 

2142.2 

-736.11 

2.12503 

.31510 

28.348 

9 

_  2 

J.B59 

1.1000 

io.r.ooo 

5000.0 

375.55 

2.23176 

.35233 

27.910 

10 

1 

.636 

1.1000 

5.O0O0 

4556.8 

140.84 

2.23175 

.35009 

28.159 

_ _ 15. 

2 

2.998 

1.1000 

l.nooo 

3421.7 

•310.77 

2.23174 

.33957 

28/545 

8 

2 

1.667 

1.1000 

.5000 

2965.6 

•465.65 

2.23174 

.33275 

28.348 

9 

2 

1.859 

1.5000 

10. 0000 

4000.0 

86.85 

2.24618 

.35716 

26.369 

4 

1 

.381 

1.5000 

s.rooo 

3462.0 

•107.61 

2.24618 

.  .  .35117 

26.376 

_  12... 

2  . 

2.432 

1.5000 

1.0000 

2450.6 

-462.64 

2.24617 

.33407 

26.377 

9 

2 

1.858 

_  1.5000 

■  5000 

2106.0 

-581.64 

_2, 2.46.03.. 

.32609 

26.377 

8 

_ 2 _ 

1.664 

1.5000 

10.0000 

5000.0 

495.39 

2.33690 

.36468 

26.231 

5 

1 

.424 

1.5000 

5 . "000 

4416.9 

248.82 

2.33690 

.36070 

26.333 

11 

2 

2.241 

1.5000 

1.0000 

3176.9 

-208.68 

2.33689 

.34721 

26,376 

11 

2 

2.241 

1.5000 

.5000 

2738.0 

-362,77 

2.33609 

,.33984 

26.377 

_  11 _ 

_ 2 _ 

2.244 

AEDC-TR-71-256 


TABLE  III 


SPECIFY  TEMPERTURE  EQUIVALENCE  RATIO  AND  PRESSURE  >o"r  A  CH2  .HYDROCARBON  WITH  AIR 
EQULI0RIUM  ISENTROPIC  EXPANSION  FROM  ID  ATM. 


.  .  _  .  EO 

p 

T 

.  SPECIES  MOL  . 

fractions 

RATIO 

ATM 

DEQ  R 

H 

0 

OH 

NO 

CO 

•  900 

10.00 

4o  0.0 

.1888-04 

.3029-03 

.2095-02 

.9194-02 

.4774-03 

.  .900 

5.00 

3485.0 

•18B6-05 

.9628-04 

.6885-03 

,4526-02  .  . 

.7285-04 

•  900 

1.00 

2463.1 

.74B9-09 

.1806-06 

.1576-04 

.4438-03 

.1190-06 

.....  .900 

.50 

21"4.4 

.6822-11 

.5727-08 

.  .1647-05 

.1144-03  . 

.2477.-08..  .  . 

.900 

10.00 

50n0 .0 

.8477-03 

.4647-02 

.1123-01 

'  .2324-01 

.8534-02 

..  .900 

5.00 

4491.9 

.2550-03 

•1912-02 

.  .>6301-02. 

_ .1535-Jll _ 

...3362-02 _  ._ 

.900 

1.00 

3319.7 

.2120-05 

.5735-0* 

.'62*0-03 

.3431-02 

.6867-04 

...  _  .900 

.50 

2865.5 

.  ..9434-07 

.SSSBtOS.. 

.1389-03. 

_ .1353^02  _ 

.  .5378-05. . .  . . 

.900 

10.00 

40'  0.0 

.3795-04 

.1322-03 

.1837-02 

.3964-02 

.1895-02 

_  ..  .900 

5.00 

35  9.6 

.4417-05 

.2739-0* 

.6462-03 

.2024-02 

.3300-03 

.900 

1.00 

2515.4 

.2829-08 

.1269-06 

.1864-04 

.2277-03 

.7903-06 

.900 

.50 

2162.2 

.3552-10 

.5083-08 

.2267-05 

,6392-04 

..2140-Q7.  ....  . 

.900 

10.00 

Sjnu.O 

.1556-02 

.2449-02 

.1085-01 

.1205-01 

.2522-01 

.900 

9.00 

4565.2 

.6407-03 

.1130-02 

.6586-02 

.7860-02 

.1390-01 

.900 

1.00 

3912.6 

.1500-04 

.6230-04 

.9751-03 

.20*1-02 

.7429-03 

.900 

.50 

3061.1 

.1026-05- 

.8757-05 

.2665-03 

,9049-03 

.8288-0* 

1.000 

10.00 

40 '0.0 

.6637-04 

.4750-0* 

.1154-02 

.1420-02 

.5647-02 

1.000 

5.00 

3536.2 

.1179-04 

.6608-05 

.3386-03 

.4510-03 

.1908-02 

1.000 

1.00 

2556.4 

.2623-07 

.5899-08 

.4437-09 

.B4B1-05 

.3857-04 

1-000 

.50 

2201 .6 

.6913-09 

.  .8092-10.. 

_  .. 3202-06.. 

.7579-06 

.3864-05 

1.000 

10.00 

5ono  .0 

.1894-02 

.1802-02 

.9731-02 

,8839-02 

.3485-01 

_  1-000 

5.00 

4577.7 

.8638-03 

.7594-03 

.5661-02 

.5154-02 

.2261-01 

1.000 

1.00 

3590.9 

.4590-04 

.2697-0* 

.7044-03 

.7057-03 

.3780-02 

1.000 

.50 

3159.5 

.5765-05 

.2526-05 

.  1623-03 

.1811-03 

.1022-02 

1.100 

10.00 

40'  0.0 

.1703-03 

.7437-05 

.4639-03 

.2199-03 

.3144-01 

.  1.100 

5.00 

3483.8 

.4176-04 

.2208-06 

.6036-04 

±1732-04 

•  2994-01. 

1.100 

1.00 

2488.1 

.4340-06 

.2341-11 

.8030-07 

.4838-08 

.2502-01 

--  i.100 

.50 

2142.2 

.2983-07 

.3056-1* 

_,l734-08 

.  .  .  ,4411-10  .. 

._lU87;01_ 

1.100 

10.00 

5C0.0 

.2384-02 

.1210-02 

.8225-02 

.5894-02 

.4957-01 

1.100 

_ 5 L.  OS . 

455.6  ,.8 

. 1146-02 

.3688-03 

.4532-02 _ *2576-02 _ 

.3911-01 

1.100 

1.00 

3421.7 

.7308-04 

.6083-06 

.9623-04 

.2571-04 

.2977-01 

.1.100 

.50 

2965.6 

....  1257.-04 

.7278-08 

_ .7435-05 

.1075-05 

.2802-01 

1.900 

10.00 

4o:o.q 

.4586-03 

.9351-06 

.1571-03 

.2634-04 

.1209*00 

.  _ 1.500 

5.00 

3462.0 

.9979-04 

.2356-07 

_  •  l8_32-_P_4  _ 

jJ878-05 

.1168*00 

1.500 

1.00 

2450.6 

.7489-06 

.1749-12 

.1873-07 

.4290-09 

.1016*00 

1.500 

.,50_ 

.210.6,0 

..4  448iOL 

.jJ9.1A-.15_ 

_  jJ5J0  J.9 

_ >3517_-U _ 

.9166-01 

1.500 

10.00 

5000.0 

.5189-02 

.2508-03 

.3711-02 

,1172-02 

.1259*00 

1.500 

5..00  . 

4416,9 

,.2026-02.. 

_ ..265.5-04.. 

>1001-02 

.,2270-03 

.1230*00 

1.500 

1.00 

3176.9 

.6433-04 

.6367-08 

,7727-05 

.5710-06 

.1139*00 

1.500 

.50 

2738.0 

.8160-05 

4354-10 

_>4?25-06 

.1655-07 

.1075*00 

AEDC-TR-71-256 


TABLE  III  (Concluded) 


SPEC! PY  TEMPER TURE  EQUIVALENCE  RaTI0'‘a’nD~ PRESSURE  FOR  A  CHE  HYDROCARBON  "il I TH  AIR" 


EQULIBRIUM  ISENTROPIC  EXPANSION  FROH  10  ATM, 

_  EQ  •  P  T  .  SPECIES _ MOL _ FRACTIONS 


RATIO 

ATM 

OEQ  R 

H2 

02 

H20 

C02 

N2 

A 

.500 

10.00 

40'  0.0 

.8844-04 

.9613-01 

.6626-01 

.6722-01 

.7492*00 

.9001*02 

_  .500 

5.00 

3485.0 

.1669-04  . 

.987.7.701 

_  ..67.10-01. 

.6768-01 

..7521*0.0 _ 

.9007-02 

.500 

1.00 

2483.1 

.5815-07 

.1010*00 

.6747-01 

.6777-01 

.7543*00 

.9009-02 

.500 

.50 

21T.4.4 

.1962-08 

,1012*0.0 

.6747-01 

.6777-01 

.7545*00 

.9009-02 

.500 

10.00 

50"0.0 

.1269-02 

,8850-01 

.5945-01 

.5851-01 

.7349*00 

.8913-02 

.500  5.00 

4491.5 

_._J542o_-.0_3  . 

■9255-01  . 

_ -6334-0.1_ 

.6400-01 

.7433*00 

.8967-02 

.500 

1.00 

3319.7 

.1711-04 

.9933-01 

.6713-01 

,6768-01 

,7526*00 

.9007-02 

.  500 

.50 

2865.5  . 

.1804-05  _ 

...,1005*00 _ 

.  ,.6740-01.. 

.6776-01 

.  ._7?38*_D0 _ 

^900.^.02  ... 

.900 

10.00 

40  10.0 

.3571-03 

.1831-01 

.1168*00 

.1165*00 

.7315*00 

.8758-02 

.900 

..5.  00.. 

.  3509,6 

.7512-04. 

.  .,1872-01. 

.  .117.8*00 

_J1»2*0D 

■7334*00 _ 

.8749-02 

.900 

1.00 

2515.4 

.3654-06 

.1959-01 

.1183*00 

.1185*00 

.7346*00 

.8773-02 

.900 

.50 

.2162.2  .. 

.  .1550-07 

__._1?68-01 

.1183*00 

.1185*00 _ 

.7347*00 

_ .8773-02 

.900 

10.00 

50,’0.0 

.4277-02 

.2450-01 

.1054*00 

.9097-01 

.7141*00 

.8599-02 

.900 

5.00 

4565.2 

.2392-02 

.2166-01 

.-.111.0*00. 

..1034*00 

.7228*00 

.8678-0? 

.900 

1.00 

3512.6 

•  1691-03 

.1886-01 

.1175*00 

,1177.00 

.7331*00 

.8766-02 

.  _  .  .900 

.50 

3061.1 

.2423-04 

,1924-01.  . 

.1161*00 

...1184*00 

.7342*00 

.8772-02 

1.000 

10.00 

40'  0.0 

.1092-02 

.2364-02 

.1284*00 

.1247*00 

.7265*00 

.8683-02 

l.noo 

.5.00 

3536.2. 

.  .4337-03 

.8603-03  .. 

.  ,.1298*00. 

.  1288*00 

. 7286*00  .8704-02 

1.000 

1.00 

2556.4 

.1711-04 

.2213-04 

.1305*00 

.1308*00 

.7299*00 

.8715-02 

1.000 

.50 

22‘<1 . 6 

.2641-05 

.2520-05 

.1305*00 

.1308*00 

.7299*00 

.8715-02 

1.000 

10.00 

SOfO.O 

.6333-02 

.1331-01 

.1151*00 

.9266-01 

.7070*00 

.8494-02 

1.000 

5.00 

4577.7 

.4122-02  . 

.9210-02 . 

.1211*00 

.J061.00 

.7158*00 

.8578-02 

1.000 

1.00 

3590.9 

.8467-03 

.1782-02 

.1290*00 

.1267*00 

.7277*00 

.8693-02" 

l.opo 

.50 

3159.5 

.2824-03 

.5209-03 

.1301*00 

.1297*00 

.7293*00 

.8709-02 

1.100 

10.00 

40:>0.0 

.7195-02 

.5796-04 

.1324*00 

.1087*00 

.7109*00 

.8489-02 

..  l.lOfl 

5.00 

3483.8 

.8259-02  ' 

.1536-05 

•1316*Q0 

.1103*00 

.7113*00 

.8493-02 

1.100 

1.00 

2488.1 

.1315-01 

.1085-10 

.1268*00 

.1152*00 

.7113*00 

.6494-02 

_ 1.100 

.50 

2142.2 

.1669-01 

.1149-13  . 

.1233*00 

.1187*00 

.7113*00 

.8494-02 

1.100 

10.00 

50r  0 . 0 

.1004-01 

.6000-02 

.1224*00 

.8849-01 

.6974*00 

.8362-02 

1.100 

5.00 

4556.8 

.8049-02 

.2430-02 

,.1284*00 

.1002*00 

,7053*00 

.8437-02 

1.100 

1.00 

3421.7 

.8440-02 

.4151-05 

.1314*00 

.1104*00 

.7113*00 

.6493-02 

1.100 

.50 

2965.6 

.1015-01 

.4233-07 

.1298*00 

.1122*00 

.7113*00 

.8494-02 

1.500 

10.00 

40*0.0 

.5217-01 

,9163-06 

.1207*00 

.5254-01 

.6454*00 

.7707-02 

1.500 

5.00 

3462.0 

.5632-01 

.2136-07 

.1168*00 

,5661-01 

.6456*00 

.7709-02 

1.500 

1.00 

2450.6 

.7158-01 

.1199-12 

.1016*00 

.7183-01 

.6456*00 

,7709-02 

1.500 

.50 

2106.0 

.8155-01 

.1108-15 

.9166-01 

j.8180-01 

.6456*00 

.7709-02 

1.500 

10.00 

50'0.0 

.4754-01 

.2578-03 

.1203*00 

.4659-01 

.6415*00 

.7666-02 

1.500 

5.00 

4416.9 

, 4967-01 

,2703-oR.. 

.  ...1217*00 

....5012-01 

•  6444  +  00 

,7696-02 

1.500 

1.00 

3176.9 

.5930*01 

.5456-08 

.1139*00 

.5958-01 

.6456*00 

.7709-02 

1.900 

.50 

2738.0 

.6567-01 

.3318-10 

.1075*00 

.6592-01 

.6456*00 

.7709-02 
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TABLE  IV 


SPECIFY  TEMPER TUrE"  EQUIVALENCE. RATIO  AND  PRESSURE  FOR  A  CH2  hydrocarbon  with  air 


FR02EN  I SENTRQP I C  EXPANSION  FROM  10  ATM. 


EQ 

p 

T_ 

H 

S 

CP 

MOL 

NO 

TYPE 

JUNE 

ratio 

ATM  DES  R 

BTU/LB 

'  RTU/LB  OR 

btu/lb  or 

HT 

ITER 

case 

SEC 

.3000 

10.0 

4000.0 

398.84 

2.06887 

.32330 

28.906 

3 

1 

.374 

.3000 

5.0 

3448...Z _ 

221.81 _ 

2.06887 

.31609  . 

28.906 

4 

3 

.463 

.5000 

1.0 

2413.4 

-99.44 

2.06888 

.30241 

26.906 

5 

3 

.343 

.3000 

■  5 

2039.2 

-2nS . 77 

_2j.a6B.89 

.29450 

26.906 

4 

3 

.464 

.3000 

10.0 

5000.0 

625.68 

2.16345 

.32903 

28.625 

3 

i 

,374 

_  .5000 

JJ 

4316.7 

602.06 

2.16346 

.32523 

26.623 

4 

3 

.463 

.3000 

1.0 

3043.4 

194.87 

2.16346 

.31264 

28.625 

5 

3 

.343 

..  .  .5000 

.5 

2603.4 

59 .31 

__2jJA347 

.30393 

28.625 

4 

3. 

,463 

.9000 

10.0 

4000. 0 

-40.11 

2.06972 

.33940 

26.662 

3 

1 

.373 

.9000 

5.0 

.  3.471,6 _ 

-217.99 

_ 2...0.8.973 

.33363 

28.862 

4 

3. 

.463 

.9000 

1.0 

2471.3 

-544.36  . 

2.06974 

.31703 

26.862 

3 

3 

.343 

.9000 

2122.0 

-651  -.7  4 

2.08973 

.30843 

28.916 

4 

3 

.464 

.9000 

10.0 

5000.0 

465.69 

2.20149 

.34521 

26.337 

4 

1 

.423 

.9000 

3.0 

4340.2 

239.20 

2.20149 

.34113 _ 

28.337 

4 

.3 

.464 

.9000 

1.0 

3099.4 

•176.86 

2.20130 

.32766 

28.337 

5 

3 

•  343. 

.  .  ..  .9000 

.3 

2668.2 

-316.72 

2...20131 

.32046 

26.337 

4 

3 

.464 

1.0000 

10.0 

4000.0 

•133.00 

2.09516 

.34316 

28.797 

5 

1 

.475 

1.0000 

3.6'  5476.0 

-311.39 

2.09317 

.33730 

28.797 

4 

3 

.463 

1.0000 

1.0 

2477.3 

-631.23 

2.09316 

.31990 

29.041 

3 

3 

.546 

...  .1..0000 

.3 

2128.9 

-728.91 

2.09520 

.31074 

29;361 

3 

.464 

1.0000 

10.0 

3000.0 

403.11 

2.21416 

.34690 

28.172 

6 

1 

.526 

.  _  ...  1.0000 

5.0 

4343,2 

177.23 

2.21418 

.34476 

26.172 

3 

.464 

1.0000 

1.0 

3106.2 

-241.86 

2.21419 

.33123 

26.172 

5 

3 

.344 

1.0000. 

-3 

2*79.7 

-374.67 

2 . 21429 

.32338 

28.386 

1 

3 

.222 

1.1000 

10.0 

4000.0 

-106.37 

2.12503 

.34596 

26.334 

3 

1 

.423 

1.1000. 

....  3.0 

. _3*7JLt0_. 

•289,38 

2.12503 

.34013 

28.334 

4 

3 

.463 

1,1000 

1.0 

2437,3 

•631.31 

2.12505 

-.32233  - 

26.340 

5 

3 

.343 

1.1000 

j.5 

.  —IMP..  3 

•766.96 

2.12510 

.31303 

26.360 

3 

3 

.345 

1.1000 

10.0 

5000.0 

373.33 

2.23176 

.35233 

27.910 

10 

1 

.636 

1 . 1D00 

3.0 

4343.4 

145.31 

2.23177 

.34606 

27.910 

4 

3 

.463 

1.1000 

1.0 

3106.3 

-277.53 

2.23178 

.33425 

27.910 

5 

3 

.343 

.U1000 

l5 

2343.3 

•456.82 

2.23178 

_ ,.32632 

27.930 

6 

3 

,626 

1.3000 

10.0 

4000.0 

86.65 

2.24616 

.35716 

26,369 

4 

1 

.361' 

1.5000 

.5.0 

...  .3431.7 

-107.31 

2.24619 

.35077 

26.369 

4 

3 

.464 

1.5000 

1.0 

2420.9 

•460.06 

2.24620 

.33165 

26,369 

5 

3 

.544 

1.3000 

.5 

2063.9 

•376.91 

2.24621 

.32241 

26.369 

4 

3 

.464 

1.3000 

10.0 

5000.0 

493.39 

2.33690 

.36466 

26.231 

i 

.423 

1.5000 

3j.Q 

4326.0 

251.11 

2.33690 

.35994 

26.231 

4 

3 

.464 

1.5000 

1.0 

3062.7 

-195,14 

2.33691 

,34474 

26.231 

5 

3 

.544 

1.5000 

.5 

2623.5 

-344.11 

2.33692 

.33642 

26.231 

4 

3 

.463 
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TABLE  IV  (Continued) 


SPECIFY  TEMPERTURE  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  A  CHI  HYOROCaRBON  WITH  AIR 


<J\ 

O 


FROZEN  I SENTROP I C  EXPANSION  FROM  10  ATM. 

EO  P _ T _ SPECIES _ HflL _ ESACHONS 


RATIO 

ATM 

DEQ  R 

H 

0 

DH 

NO 

CO 

•  300 

10-00 

4000.0 

■1888-04 

.3029-03 

.2095-02 

.9194-02 

.4774-03 

.Snn 

5.00 

3440.2 

.1888-04 

.3029-03 

.2095-02 

.9194-02 

_ ,.4774-03 

•  500 

1.00 

2415.4 

.l8SS'04 

.3029-03 

.2095-02 

.9194-02 

.4774-03 

.500 

.50 

2059.2 

,1688-04. 

.3079-03 

.2095-02 

.9194-02 

.4774-03 

.500 

10.00 

SOOO.O 

.6477-03 

.4647-02 

.1123-01 

.2324-01 

.8534-02 

.9nn 

5.00 

4316.7 

.6477-03 

. 4647-02 

.1123-01 

.2324-01 

.8534-02 

.500 

1.00 

3043.4 

.8477-03 

.4647-02 

. .1123-01 

.2324-01 

.6534-02 

.500 

.50 

2605  .A. 

.8477-03 

. 4647-02 

.1123-01 

.2324-01 

.8534-02 

•  900 

10.00 

4000.0 

.3795-04 

.1922-03 

.1837-02 

.3964-02 

.1895-02 

.9nn 

5.00 

3471  .0 

37  95-04 

.1322=03- 

.1837-02 

.3964-01 

.1895-02 

•  900 

1.00 

2471.3 

.3795-04 

-  .1322-03 

.1637-02 

.3964-02 

.1895-02 

.900 

.50 

2122.0 

.3795-04 

.1322-03. 

.1837-02 

.3964-02 

.1895-02 

.900 

10.00 

5000-0 

.1556-02 

.2445-02 

.1085-01 

.1205-01 

.2522-01 

.900 

5.00 

4340.2 

.1556-02 

■2445-02 

.1085-01 

.1209*01 

.2922=01 

.900 

1.00 

3099.4 

.1556-02 

.2445-02 

.1085-01 

.1205-01 

.2522-01 

.900 

.50 

2666.2 

.1556-02 

.2445-02 

.1085-01 

.1205-01 

.2522-01- 

1.000 

10.00 

4000.0 

.6637-04 

.4750-04 

.1154-02 

.1420-02 

.5647-02 

_ l.noo 

5.00 

3476.0 

-.6637-04 _ 

.47 an -a* 

.1154-07 

_ .14ZQ-BZ 

.3847-02 _ 

1.000 

1.00  2477.3 

.6637-04 

,4750-04 

.1154-02 

.1420-02 

.5647-02 

_ 1.9.0.0  0 

...  .50  _2128_._S_ 

.6637-04 

.4750-04 

-.1154-02 

.1420-02 

.5647-02 

1.000 

10.00  5000.0 

.1894-02 

.1802-02 

.9731-02 

.6839-02 

.3485-01 

l-lOOO. 

5.00  43_43.2. 

.1694-02 

_ i-16Q2~02 

.9731-02 

.6839-02 

.3485-01 

1.000 

1.00  3106.2 

.1894-02 

.1802-02 

.9731-02 

.8839-02 

.3465-01 

1.0  OQ 

.50  2679.7 

.1694-02 

.1802-02 

.9731-02 

.8839-02 

.3485-01 

1.100 

10.00  4000.0 

.1703-03 

.7437-o! 

.4639-03 

.2199-03 

.3144-01 

.  .  .  .1-100 

5.00..  3472. 

..  ,1703-03 

. 7437-Q5 

.  .4639-03 

.2199-03 

■  3144-01 _ 

1.100 

1.00  2437.3 

.1703-03 

.7437-05 

.4639-03 

.2199-03 

.3144-01 

_ J..100 

.50  1990.3 

.1703-03 

,7437-05 

.4639-03 

.2199-03 

.3144-01 

1.100 

10.00  5000.0 

.2384-02 

.1210-02 

.8225-02 

.5694-02 

.4957-01 

1.100 

5.00  4343.4 

.2384.02 

.1210-02 

.8225-02 

.5894-02 

.4957-01 

1.100 

1.00  3106.3 

.2384-02 

,1210-02 

.8225-02 

.5694-02 

.4957-01 

1.100 

■50  2545.3 

.2384-02 

■1210-02 

.8225-02 

,5894-92 

.4957-01 

1.500 

10.00  4000.0 

.4586-03 

.9351-06 

.1571-03 

.2634-04 

.1209*00 

1  •  soo__ 

5.00  3451.7 

.4586-03 

.9351-06 

.1571-03 

.2634-04 

.1205*00  . 

1.500 

1.00  2420.9 

.4586-03 

.9351-08 

.1571-03 

.2634-04 

.1209*00 

l.snn 

.50  2063.9 

•4986-03 

.9351-06 

.1571-03 

.2634-04 

.1209*00 

1.500 

10.00  5000.0 

.5189-02 

.2508-03 

.3711-02 

.1172-02 

•1259*00 

1.500 

5jJOO  4326. J). 

.5189-02 

.2506-03 

.3711-02 

.1172-02 

.1259*00 

1.500 

1.00  3062.7 

.5189-02 

.2508-03 

.3711-02 

.1172-02 

.1259*00 

_ 1_*_5  OQ 

.  ,?Q_?425JL 

_ .9189-02 

.2508-03 

.3711-02 

■1172-02 

.1259*00 
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TABLE  IV  (Concluded) 


SPECIFY  TEMPERTURE  EQUIVALENCE  RATIO  AND  PRESSURE  FOR  A  CHZ  HYDROCARBON  WITH  AIR 
frozen  ISEntropic  expansion  fr"oh"T6  "atm".’  ~~ 


EO 

P 

T 

SPECIES  MOL. 

_e.Ractio.ns 

RATIO 

atm 

OEQ  R 

H2 

02 

H20 

COZ 

N2 

k 

.300 

10.00 

40n0 • 0 

.8844-04 

.9613-01 

.6626-01 

.6722-01 

.7492*00 

.9001-02 

_ _ ..500 

5.00 

3448.2 

.8844-04 

.9613-01 

,  6626-01... 

_ ,6722-01  . 

7492*00 .  . 

.9001-02 

.300 

1.00 

2415.4 

.8844-04 

.9613-01 

.6626-01 

.6722-01 

.7492*00 

.9001-02 

.500 

_  .30 

2059.2 

.8844-04 

.9613-01 

...  ..6626-01... 

.6722-01 

.7492*00 

.9001-02 

.300 

10.00 

50^0.0 

.1269-02 

.8830-01 

.5945-01 

.5851-01 

.7349*00 

.8913-02 

.500 

i.QO 

4310  *7 

.1269-0? 

.BBSQ-01. 

.5943-01 

■5831-01 

.7349*00 

.8913-02 

.300 

1.00 

3043.4 

.1269-02 

.8830-01 

.3943-01 

.3831-01 

.7349*00 

.8913-02 

.300 

.  ..  .30 

26rS.4 

.1269-02 

.8850-01 

_..394J-01.._ 

_•  3831 -01 _ 

.7349*0.0 

.8913-02 

.900 

10.00 

4oro  .o 

.3371-03 

•1B31-01 

.1168*00 

.1163*00 

.7313*00 

.8738-02 

.900 

5.00 

3471.8 

..‘3371-03 

...  .1831-01 

___Ul6a.*00.... 

_ *1 165*  0.0.  _ 

■  7315*0.0  . 

.8738-02 

.900 

1.00 

2471.3 

.3371-03 

.1831-01 

.1168*00 

.1163*00 

.7313*00 

.8738-02 

.900 

.  ,5o_ 

2122.0 

.3371-03 

_i200Q-0l _  .1166-00 

.1165*00 

■7315*00 

.8758-02 

.900 

10.00 

3000.0 

.4277-02 

.2450-01 

.1034*00 

.9097-01 

.7141*00 

.8399-02 

.900 

3.00 

4340.2 

.4277-02 

.2450-01 

_ ,1034*00  _ 

...9097rj)l. 

.  ,7141*00..  _ 

.8599-02 

.900 

1.00 

3099.4 

.4277-02 

.2430-01 

.1054*00 

.9097-01 

.7141*00 

.8599-02 

.900 

.30 

2668.2. 

_. 4277-02 

...2430-OJ 

j. 1.034*00 

j9097-01 

...  7141*00  ... 

.8599-02 

1.000 

10.00 

4000.0 

.1092-02 

.2364-02 

.1284*00 

.1247*00 

. 7265*00 

. BOB  3-02 

1.000 

J? .  go 

3476.0 

.1092-02 

.2364-02 

.1284*00 

_ J247.00 

. 7263*00 

.8683-02 

1.000 

1.00 

2477.3 

.1092-02 

.1000-01 

.1284*00 

.1247*00 

.7263*00 

. B683-02 

1.000 

.50 

2128*9 

.1092-02 

.  2000-.01 

.1281*00 

.1247*00 

.7265*00 

_,66BJ-02  . 

1.000 

10.00 

50*0.0 

.6333-02 

.1331-01 

.1151*00 

.9266-01 

.7070*00 

.8494-02 

1.000 

3.00.. 

.4343.2 

.6333-02 

.  .  .1331-01. 

.-*113.1*0.0 

_ .92««-01 

_,.7Q7_P*00  . 

.8494-02 

1.000 

1.00 

31C6.2 

.6333-02 

.1331-01 

.1151*00 

.9266-01 

.7070*00 

.8494-02 

1.000 

..  __i3.0  ._ 

2679.  .7. 

.6333-02 

,.2000-01 

.1131*00 

.9266-01 

.7070*00 

.8494-02 

1.100 

10.00 

40*0.0 

.7195-02 

.5796-04 

.1324*00 

.1087*00 

.7109*00 

.8489-02 

1_.100 

3.00 

3472.0 

.7193-02 

.3796-04 

.1324*00 

.1087.00 

.7109*00 

.8489-02 

1.100 

1.00 

2437.3 

.1000-01 

.5796-04 

.1324*00 

.1087*00 

.7109*00 

.8489-02 

1.100 

.  .50 

1990.3 

.2000-01 

•  5796^0.4 

.1324*00 

.1087*00 

.7109*00 

.B469-02 

1.100 

10.00 

30~0.0 

.1004-01 

.6000-02 

.1224*00 

.8849-01 

.6974*00 

.6362-02 

1.100 

_ 3.00 

4343.4 

_ .1004-01 

.6000-02 

.1224*00 

.8849-01 

.6974*00 

■8362-02 

1.100 

1.00 

3106.3 

.1004-01 

.6000-02 

.1224*00 

.6849-01 

.6974*00 

.8362-02 

1.100 

.30 

2543.3 

.2000-01 

.6000-02 

.1224*00 

.8649-01 

.6974*00 

.8362-02 

1.500 

10.00 

40*0.0 

.5217-01 

.9163-06 

,1207*00 

.3234-01 

.6454*00 

.7707-02 

1.500 

5.00 

3451.7 

.5217-01 

.9163-06 

.1207*00 

,3254-01 

.6454*00 

.7707-02 

1.500 

1.00 

2420.9 

.3217-01 

.9163-06 

.1207*00 

.5234-01 

.6454*00 

.7707-02 

1.300 

j.59 

206.3.9 

_»??1.7-B1_ 

.9163-06 

.1207*00 

.5254-01 

.6434*00 

.7707-02 

1.500 

10.00 

SOi'O.O 

.4754-01 

.2378-03 

.1203*00 

.4639-01 

.6415*00 

.7666-02 

_ 1.300. 

3 . 00._ 

4326.0 

,4754-01 

_.  j2S7V.03 

_..  ,1203*00 

,4659- 01_ 

..,51 15 *00 

.7666-02 

•  1.300 

1.00 

3062.7 

.4734-01 

.2578-03 

,1203*00 

.4659-01 

.6415*00 

.7666-02 

1.500 

.50 

2625.5 

.4754-01 

.2578-03 

.1203*00 _ 

.6415*00 

.7666-02 

AEDC-TR-71  -256 


UNCLASSIFIED 

Security  Classification 


DOCUMENT  CONTROL  DATA  •  R  &  D 

(Security  elms  s  ificstion  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  la  clmaeltledl 


3  REPORT  TITLE 

AN  EFFICIENT  NUMERICAL  METHOD  FOR  THE  CALCULATION  OF  CHEMICAL 
EQUILIBRIUM  IN  THE  H/C/O/N/A  SYSTEM 

4  DESCRIPTIVE  NOTES  (Type  ot  report  and  inclusive  dates) 

Final  Report  -  September  to  October  1970 _ 

fl  au  THO R IS l  (First  name,  middle  Initial,  taat  name) 

I.  T.  Osgerby  and  R.  P.  Rhodes,  ARO,  Inc. 


10  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT 

An  efficient  method  is  presented  for  calculating  equilibrium  chemi¬ 
cal  composition  at  a  given  pressure  and  stoichiometry  with  temperature, 
enthalpy,  or  entropy  specified  as  the  additional  condition.  A  sub¬ 
stantial  reduction  in  computational  time  was  achieved  over  conventional 
methods  by  employing  a  linearization  technique  and  reducing  the  number 
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tions.  Examples  are  given  for  the  H2/02  system  (six  equations  reduced 
to  two)  and  C/H2/02/N2/Ar  system  (eleven  equations  reduced  to  three) . 

A  computer  program  is  presented  for  solving  problems  with  fuels  which 
involve  C  and  H  (hydrocarbon  or  hydrogen)  or  N  and  H  (such  as  hydrazine) 
and  with  oxidants  such  as  02,  H202,  and  HN03 .  Extension  of  the  method 
to  include  other  systems  is  explained  in  sufficient  detail  to  eliminate 
the  requirement  for  an  extensive  knowledge  of  chemistry  by  potential 
users. 
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